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ABSTRACT 

X-ray absorption line spectroscopy has recently shown evidence for previously unknown Ultra-fast 
Outflows (UFOs) in radio-quiet AGNs. These have been detected essentially through blueshiftcd 
Fe XXV/XXVI K-shell transitions. In the previous paper of this series we defined UFOs as those 
highly ionized absorbers with an outflow velocity higher than 10,000 km/s and assessed the statistical 
significance of the associated blueshifted absorption lines in a large sample of 42 local radio-quiet 
AGNs observed with XMM-Newton. The present paper is an extension of that work. First, we 
report a detailed curve of growth analysis of the main Fe XXV/XXVI transitions in photo-ionized 
plasmas. Then, we estimate an average SED for the sample sources and directly model the Fe K 
absorbers in the XMM-Newton spectra with the detailed Xstar photo-ionization code. We confirm 
that the frequency of sources in the radio-quiet sample showing UFOs is >35% and that the majority 
of the Fe K absorbers are indeed associated with UFOs. The outflow velocity distribution spans from 
-10,000 km/s (-0.03c) up to -100,000 km/s (-0.3c), with a peak and mean value of -42,000 km/s 
(—0.14c). The ionization parameter is very high and in the range log£— 3-6 erg s _1 cm, with a mean 
value of log£— 4.2 erg s _1 cm. The associated column densities are also large, in the range Njj~ 10 22 - 
10 24 cm" 2 , with a mean value of 7V#— 10 23 cm~ 2 . We discuss and estimate how selection effects, 
such as those related to the limited instrumental sensitivity at energies above 7 keV, may hamper the 
detection of even higher velocities and higher ionization absorbers. We argue that, overall, these results 
point to the presence of extremely ionized and possibly almost Compton thick outflowing material in 
the innermost regions of AGNs. This also suggests that UFOs may potentially play a significant role 
in the expected cosmological feedback from AGNs and their study can provide important clues on the 
connection between accretion disks, winds and jets. 

Subject headings: black hole physics — galaxies: active — galaxies: Seyfert — line: identification — 
plasmas — X-rays: galaxies 



1. INTRODUCTION 

Blueshifted Fe K absorption lines have been detected in 
recent years between 7 and 10 keV in the X-ray spectra 
of several radio-quiet AGNs (e.g., Chartas et al. 2002; 
Chartas et al. 2003; Pounds et al. 2003; Dadina et 
al. 2005; Markowitz et al. 2006; Braito et al. 2007; Cappi 
et al. 2009; Reeves et al. 2009). These lines are com- 
monly interpreted as due to resonant absorption from Fe 
XXV-XXVI associated with a zone of circumnuclear gas 
photo-ionized by the central X-ray source, with ioniza- 
tion parameter log£— 3-5 erg s _1 cm and column den- 
sity Nh~ 10 22 -10 24 cm~ 2 . The energies of these absorp- 
tion lines are systematically blueshifted and the corre- 
sponding velocities can reach up to mildly relativistic 
values of —0.2-0. 4c. These findings are important be- 
cause they suggest the presence of previously unknown 
massive and highly ionized absorbing material outflow- 
ing from their nuclei, possibly connected with accretion 
disk winds/ejecta (e.g., King & Pounds 2003; Proga & 
Kallman 2004; Sim et al. 2008; Ohsuga et al. 2009; King 
2010a; Sim et al. 2010) or the base of a possible weak jet 
(e.g., Ghisellini et al. 2004). 

In particular, a uniform and systematic search for 



blueshifted Fe K absorption lines in a large sample of 
42 radio-quiet AGNs observed with XMM-Newton has 
been performed in Tombesi et al. (2010a, hereafter TA). 
This study allowed us to assess their global detection 
significance and to overcome any possible publication 
bias. For clarity, there we defined Ultra-fast Outflows 
(UFOs) those highly ionized absorbers detected essen- 
tially through Fe XXV and Fe XXVI K-shell absorption 
lines with blueshifted velocities v>10,000 km/s (—0.03c). 
This was done to distinguish from the classical warm ab- 
sorbers (WA) in AGNs, which are commonly detected in 
the soft X-ray band of more than half of radio-quiet ob- 
jects. In fact, these absorbers are usually less ionized, are 
outflowing with velocities lower than —1000 km/s and are 
supposed to have a different physical origin (e.g., Blustin 
et al. 2005; McKernan et al. 2007). 

In TA we detected UFOs in >35% of the sources, with 
outflow velocities even larger than 0.1c in —25% of them, 
and observed also short time-scale variability of —days 
in a few cases. However, in TA, the narrow absorption 
lines were modeled with simple inverted Gaussians and 
the blueshifted velocities were estimated assuming a line 
identification essentially as Fe XXV/XXVI ls-2p and Is- 



2 



Tombesi et al. 



3p transitions. Here we extend that work and model 
the Fe K absorbers using the detailed photo-ionization 
code Xstar v. 2.2 (Kalhnan & Bautista 2001), specif- 
ically developed for X-ray astronomy, which allows to 
derive a physically self-consistent estimate of the ioniza- 
tion parameter, column density and outflow velocity of 
the UFOs. 

In addition, Tombesi et al. (2010b, hereafter TB) also 
detected, for the first time, UFOs in three out of five 
radio-loud AGNs observed by Suzaku. These are in 
the form of a series of narrow absorption lines at en- 
ergies greater than 7 keV. Their likely interpretation as 
blueshifted Fe XXV/XXVI K-shell resonance lines im- 
plies an origin from highly ionized gas outflowing with 
mildly relativistic velocities, in the range u~0. 04-0. 15c. 
A fit with specific photo-ionization models give ioniza- 
tion parameters in the range log£~4-5.6 erg s _1 cm 
and column densities of 7V#~10 22 -10 23 cm~ 2 . Their 
estimated location within ^0.01-0.3pc from the central 
super-massive black hole (SMBH) suggests a likely origin 
related to accretion disk winds/ejecta and maybe the jet 
present in these sources. 

Depending on the absorber covering factor and geom- 
etry, the mass outflow rate of these UFOs can be com- 
parable to the accretion rate and their kinetic power can 
correspond to a significant fraction of the bolometric lu- 
minosity and is comparable to the typical jet power of 
radio-loud sources. Therefore, they might have the pos- 
sibility to bring outward a significant amount of mass and 
energy, which can have an important influence on the sur- 
rounding environment (e.g., see review by Cappi 2006). 
In fact, feedback from the AGN is expected to have a 
significant role in the evolution of the host galaxy, such 
as the enrichment of the ISM or the reduction of star for- 
mation (e.g., see review by Elvis 2006 and Fabian 2009). 
Moreover, the ejection of a substantial amount of mass 
from the central regions of AGNs can also inhibit the 
growth of the SMBHs and could also explain some fun- 
damental relations (e.g., King 2010a, b, and references 
therein), potentially affecting their evolution. Therefore, 
these UFOs can play a significant role in the expected 
cosmological feedback from AGNs and can give us further 
clues on the relation between the accretion disk and the 
formation of winds/jets in both radio-quiet and radio- 
loud sources. 

This paper is structured as follows. In §2 we describe 
a curve of growth analysis of the main Fe XXV/XXVI 
K-shell absorption lines. In §3 we present the calculation 
of the average SED of the TA sample sources and the fits 
to the XMM-Newton spectra with the photo-ionization 
code Xstar. In §4 we discuss the results of the fits and 
show the distribution of the ionization parameter, col- 
umn density and outflow velocity of the UFOs. The con- 
clusions are reported in §5. In Appendix A we report 
the detailed photo-ionization modeling for each source 
and in Appendix B the x 2 statistic and absorber redshift 
contour plots for each observation. We refer to a succes- 
sive paper III of this series (Tombesi et al. in prep.) for 
a detailed discussion of the physical properties, geome- 
try and energetics of the UFOs here analyzed and also a 
comparison with the classical warm absorbers. 

2. CURVE OF GROWTH ANALYSIS OF THE FE K 
ABSORPTION LINES 



Due to the intense X-ray luminosity of the central 
source, absorbers in the inner regions of AGNs are 
mainly photo-ionized by this radiation instead of being 
locally ionized by thermal collisions (e.g., Kinkhabwala 
et al. 2002; Crenshaw et al. 2003). In this way, they can 
reach very high levels of ionization maintaining relatively 
low and constant temperatures of T^10 5 -10 6 K (e.g., 
Nicastro et al. 1999; Bianchi et al. 2005). In order to 
properly parameterize the observed absorption features 
we need first to have a clear quantitative idea of the gen- 
eral dependence and effects of the different physical pa- 
rameters of the photo-ionized plasma on the strength of 
the Fe K absorption lines. These effects and dependen- 
cies have already been studied by Bianchi et al. (2005) 
and Risaliti et al. (2005). Nevertheless, these authors 
were mainly interested in low velocity photo-ionized out- 
flows in AGNs, instead of the ultra- fast outflows consid- 
ered here, and used slightly different assumptions. How- 
ever, the radiation processes are essentially the same in 
both cases and their results can be directly compared 
with ours. Therefore, in this section we describe a de- 
tailed curve of growth analysis of the highly ionized Fe 
K-shell absorption lines and then compare with what ob- 
served in the XMM-Newton spectra of TA. 

2.1. Fe XXV and Fe XXVI K-shell lines 

We consider the most intense K-shell transitions of Fe 
XXV and Fe XXVI reported in Table 1. All line param- 
eters throughout this paper have been taken from the 
NISlQ atomic database, if not otherwise stated. 

The most prominent Fe XXV lines are due to the Hea 
(ls 2 -ls2p) and He/3 (ls 2 -ls3p) transitions, which further 
subdivide into inter-combination (i) and resonance lines 
(r). At the moderate energy resolution of the XMM- 
Newton European Photon Imaging Camera (EPIC) pn 
instrument (FWHM>150 eV at 6 keV) these fine struc- 
ture line components are not distinguishable and they are 
blended. Therefore, we averaged their energies, weighted 
for the respective oscillator strengths. The consequent 
line energies are 6.697 keV for the Fe XXV Hea and 
7.880 keV for the Fe XXV He/3, respectively, see Table 1. 
The most intense lines from Fe XXVI are instead due 
to the Lyman series, in particular the Lya (ls-2p) and 
Ly/3 (ls-3p), each of which further subdivide into two 
resonance doublets, ri and T2- Their average energies, 
weighted for their oscillator strengths, are 6.966 keV for 
the Fe XXVI Lya and 8.250 keV for the Fe XXVI Ly/3, 
respectively, see Table 1. 

It should be noted that these strong lines are always 
accompanied by further less intense, higher order K- 
shell resonances, such as: Fe XXV He7 (ls 2 -ls4p) at 
E-8.295 keV, Fe XXV ReS (ls 2 -ls5p) at E-8.487 keV, 
Fe XXVI Ly 7 (ls-4p) at E-8.701 keV and Fe XXVI 
Ly<5 (ls-5p) at E~8.909 keV. However, their detection is 
hampered by the fact that they are intrinsically weaker 
than the former, with oscillator strengths <5% and there- 
fore with less intensity, and they are located at E>8 keV 
where the S/N and energy resolution of the EPIC pn 
spectra are worse. Therefore, here we do not present a 
curve of growth analysis of these weak lines but they are 
consistently taken into account when fitting the spectra 
using the Xstar code (see §3.2). 

1 http:/ /physics. nist.gov/PhysRefData/ASD/index. html. 
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2.2. Theoretical calculation 

We performed a curve of growth analysis of the most 
intense Fe XXV/XXVI K-shell absorption lines gener- 
ated in a photo-ionized plasma. We used the line param- 
eters reported in Table 1. This study allows us to derive a 
quantitative estimate of the dependence of the equivalent 
width (EW) of the lines, directly observable in the X-ray 
spectra, on the main physical parameters describing the 
absorbing material, such as the ionization parameter £, 
the column density Nh and the velocity broadening pa- 
rameter b. All calculations were performed separately for 
each line component listed in Table 1 and only their final 
EWs were summed. 

The relative abundances of Fe XXV and Fe XXVI ions 
are calculated by means of the photo-ionization code Xs- 
tar as a function of the ionization parameter and col- 
umn density. We assume a power-law ionizing continuum 
with r=2, which is the typical value for Seyfert galax- 
ies in the X-rays (e.g., Bianchi et al. 2005; Dadina 2008 
and references therein). We use the standard definition 
of the ionization parameter £ = Li on /nr 2 , in units of 
erg s _1 cm (Tarter, Tucker & Salpeter 1969), where Li on 
is the ionizing luminosity between 1 Ryd and 1000 Ryd 
(1 Ryd=13.6 eV), n is the number density of the ma- 
terial and r is the distance of the gas from the central 
source. This expression directly derives from the ion- 
ization balance equations and essentially represents the 
ratio between the density of the ionizing photons and the 
density of the absorbing material. 

We assume standard Solar abundances and in partic- 
ular an iron abundance relative to hydrogen of Ap c = 
3.16 x 10~ 5 (Asplund et al. 2009). We assume a constant 
density of the gas, with the typical value used for the 
Broad Line Region of type 1 AGNs, n—10 10 cm~ 3 (e.g., 
Crenshaw et al. 2003 and references therein). We did 
not fix the gas temperature but let the code calculate 
the thermal equilibrium. However, the resultant tem- 
peratures are always consistent with the typical value of 
T-10 6 K (e.g., Nicastro et al. 1999; Bianchi et al. 2005). 
Finally, we limited the treatment to a Compton thin ab- 
sorber, Nh<10 24 cm~ 2 , distributed on a geometrically 
thin shell, Ar/r = A^ 1 ^.^)-!^ < 1. 

These assumptions are slightly different from those of 
Bianchi et al. (2005) and Risaliti et al. (2005). In partic- 
ular, they derived the iron ions populations using a dif- 
ferent photo-ionization code, CLOUDY (Fcrland 2003), 
but considered the same power-law continuum with T=2 
as in this work. They assumed a lower gas density of 
n=10 6 cm~ 3 and fixed the gas temperature to T=10 6 K. 
They considered a slightly different iron abundance rela- 
tive to hydrogen of Ap c = 4.68 x 10~ 5 (Anders & Grevesse 
1989). They also used a different definition of the ioniza- 
tion parameter, U Xl which however can be simply related 
to ours as log£~4+log£/ x , for the T=2 power-law contin- 
uum considered here. 

To derive the EW of the absorption lines we followed 
an approach similar to that of Nicastro et al. (1999) and 
Bianchi et al. (2005). The EW of a resonant absorption 
line is defined as: 

r +oo 

EW= / (l-e- T ")d^, (1) 
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where t v is the dimensionless frequency specific optical 
depth of the considered transition: 



t v = / a,As = niL fiu®{v) • (2) 

Jo Tn e c 

In the above equation a v is the absorption coefficient 
at the frequency v, n>i is the number density of ions of the 
given element that populate the lower level of the transi- 
tion considered, L is the linear extension of the material 
along the line of sight, /;„ is the oscillator strength of 
the electron transition from the lower to the upper level 
and $(v) is the normalized Voigt profile of the line. We 
then directly integrated the Voigt profile, which can be 
written as (e.g., Rybicki & Lightman 1979): 

m = ^H M , (3) 
where H(a, u) is the Voigt function: 

a f + °° e~ y2 
H(a,u) = - - T — -dy. (4) 

W-oo a 2 + (u-y) 2 

in which a — ^/(AttAud) and u = (y — vq)/ Avr>. 

This is essentially a convolution of the Gaussian and 
Lorentzian line profiles. The Gaussian profile, which rep- 
resents the core of the line, depends on the Dopplcr 
broadening Avd = fob/c, where vq is the rest-frame 
specific frequency of the transition and b is the veloc- 
ity broadening parameter. The latter is the square root 
of the quadratic sum of the thermal velocity of the ions, 
v 2 h = 2kT/m, where T is the gas temperature and m 
the ion mass, and the gas turbulent velocity, or velocity 

1 /2 

dispersion along the line of sight, v tu : b — (v 2 h + v 2 u ) 

The expected temperature of a highly ionized gas in 
photo-ionization equilibrium in the central regions of 
AGNs is of the order of T-10 6 K (e.g., Bianchi et al. 2005 
and references therein) and the opacity is dominated by 
Fe ions. We can estimate that the consequent thermal 
velocity of the gas is less than 20 km/s and therefore 
its contribution to the velocity broadening is negligible, 
less than 2%, for a typical value of the turbulent veloc- 
ity of ^1000 km/s (e.g., Nicastro et al. 1999; Bianchi et 
al. 2005). Therefore, the velocity broadening parameter 
is dominated by the turbulence component and this fact 
is also strengthened by the linear dependence b oc v tu , 
instead of 6 oc T 1 / 2 . On the other hand, the Lorentzian 
profile contributes to the broadening of the line wings 
and is a function of the natural broadening parameter 
7, which in turn depends on the Einstein coefficient for 
spontaneous emission for the transition considered, A u i 
(e.g., Rybicki & Lightman 1979). 

The curve of growth of the Fe XXV Keen and 
Fe XXVI Lya lines for a fixed ionization parameter 
log£=3.5 erg s _1 cm are reported in the upper and lower 
panels of Fig. 1 , respectively. For relatively low column 
densities, less than ~10 22 cm~ 2 , the curves follow a lin- 
ear relation and the increase in EW is proportional to 
Nh, EW oc t oc Nnfiu- For larger columns, instead, 
saturation effects become more intense and the increase 
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in EW is lower. However, here we can clearly see the ef- 
fects of increasing the velocity broadening, from 20 km/s 
up to 10,000 km/s. The Gaussian line core broadens and 
the column density needed to saturate the line is larger. 
Therefore, the linear branch of the curve of growth spans 
over a larger range of Nh before saturating and the lines 
can reach higher EWs. In fact, we can see that for col- 
umn densities of V H ~10 23 -10 24 cmT 2 the EW of both 
lines can reach values of ~10 eV for 6~20 km/s and up 
to a few -100 eV for 6~10,000 km/s. 

As already noted before, given that the temperature 
of the photo-ionized gas remains almost constant at the 
value of T-10 6 K (e.g., Nicastro et al. 1999; Bianchi et 
al. 2005), which corresponds to a thermal velocity of the 
Fe ions lower than ~20 km/s, the broadening is essen- 
tially dominated by an increase in the turbulence com- 
ponent. We note that absorption lines with turbulent 
velocities as high as ~1000 km/s are shown by many 
broad-line AGNs (e.g., Crenshaw et al. 2003 and ref- 
erences therein) and such values are easily reached in 
gas moving close to SMBHs (e.g. Nicastro et al. 1999; 
Bianchi et al. 2005 and references therein). 

The dependence of the EW of the Fe XXV Hea and 
Fe XXVI Lya absorption lines with respect to the ion- 
ization parameter is shown in Fig. 2 for a set of given 
column densities and velocity broadening. The location 
of the peak of the distributions depends on the specific 
ionic population for different ionization states. In fact, 
the highest abundance for Fe XXV is reached around 
log£~3 erg s" 1 cm. Consequently, Fe XXVI peaks at 
higher ionization levels of log£~3.5-4 erg s _1 cm. More- 
over, it should be noted that the EW distribution of 
Fe XXV Hea is narrow, with a sharp rise and de- 
cline between log£~2.5-4 erg s _1 cm. Instead, that 
of Fe XXVI Lya is much broader, with a sharp rise 
at log£~3 erg s~ 1 cm and a tail extending even up to 
log£~6 erg s" 1 cm. Therefore, the Fe XXVI transi- 
tions clearly dominate for ionization levels higher than 
log£~4 erg s" 1 cm. It is also important to note the rise 
of the EW following an increase of the velocity parame- 
ter, which is equivalent to say for higher turbulent veloc- 
ities. In fact, increasing the velocity broadening b from 
100 km/s up to 5000 km/s, the resultant value at the 
peak can reach values from a few ~10 eV up to a few 
-100 eV for V H -10 23 -10 24 cm" 2 (see Fig. 2). This, 
again, shows that the increase in turbulent velocity is an 
efficient mechanism to reduce the line saturation. 

Fig. 3 shows the ratio between the EWs of the Fe 
XXV Hea and Fe XXVI Lya absorption lines as a func- 
tion of the ionization parameter and for different column 
densities. We can note that for an ionization level of 
log£~3.5 erg s _1 cm the ratio is ~1 and therefore we 
expect to observe both lines in the spectrum with com- 
parable EWs. Instead, for lower or higher values of the 
ionization parameter the Fe XXV Hea or the Fe XXVI 
Lya respectively dominate. In these cases we expect to 
observe essentially only one of these two strong lines. 
However, we can note that increasing the column den- 
sity the ratio tends to approach unity for a wider range of 
ionization levels. In the limit of Vff^lO 24 cm -2 , there is 
a plateau in the range log£~3-4 erg s _1 cm in which the 
ratio is ~1. Therefore, also in this case we would expect 
to observe both lines in the X-ray spectrum with similar 
EWs. This is again an effect due to the line saturation 



and it is less intense for higher velocity parameters. 

The EW ratios of Fe XXV Hea/He^ and Fe XXVI 
Lya/Ly/3 for different column densities and velocity pa- 
rameters are shown in the upper and lower panels of 
Fig. 4, respectively. The ionization level is fixed to 
log£=3.5 erg s _1 cm for reference. We can note that 
for column densities lower than ~10 21 -10 22 cm -2 the ra- 
tios are essentially dictated by the ratios of the different 
oscillator strengths of the transitions, ^5 (see Table 1). 
However, there is a tendency to approach unity for higher 
column densities. This effect is again due to the line sat- 
uration. In fact, for a fixed velocity broadening, the tran- 
sitions with higher oscillator strengths, in this case Hea 
and Lya, quickly reach high EWs and therefore their 
cores saturate earlier than for the weaker lines, whose 
EWs can instead still grow. It follows that for high col- 
umn densities, Vff>10 22 cm -2 , there is the possibility 
to observe both Hea/Hc/3 and Lya/Ly/3 absorption lines 
with comparable EWs in the spectrum. This effect is 
however mitigated by an increase in the velocity broad- 
ening. As clearly shown in the upper and lower panels 
of Fig. 4, higher values of b cause the decrease trend of 
the ratios to shift toward higher columns and the initial 
values are maintained for a wider range of Nh- 

Finally, we investigate the effects of changing differ- 
ent input parameters to the Xstar code on the curve of 
growth of the Fe K absorption lines. First, we checked 
the effects of assuming different continuum slopes. We 
find that changes of the power-law photon index in the 
wide range of values r=1.5-2.5 provide EWs which are 
always consistent within 10%. In particular, this es- 
sentially causes a slight shift of the peak of the EW 
distribution towards higher/lower ionization parameters 
for flat /steep slopes. This effect is mainly due to the 
fact that, given the same bolometric luminosity, in the 
flat /steep power-law cases we have more/less high energy 
photons at E>6 keV that can effectively ionize iron atoms 
and produce Fe XXV /XXVI absorption lines. There- 
fore, to obtain the same iron ionic populations, steep 
spectra require slightly higher ionization levels than flat 
ones. From the definition of the ionization parameter, 
and assuming constant luminosity and density, this can 
be translated into the need for the material to be closer 
to the X-ray radiation source for steep spectra in order to 
get more hard X-ray flux to increase the Fe XXV/XXVI 
populations. 

We also tested that the inclusion of a possible UV- 
bump between ~ 1-100 eV in the SED yields negligible 
effects to the Fe K band based results. In fact, as al- 
ready stated by McKernan et al. (2003), the presence or 
absence of the bump in the SED does not change sig- 
nificantly the parameters of the photo-ionized gas in the 
Fe K band because the main driver here is again the 
ionizing continuum in the X-rays. We checked that the 
gas density does not play an important role in the over- 
all modeling and the results are basically the same for a 
wide range of values, n~10 2 -10 12 cm -3 (e.g., Yaqoob et 
al. 2003; Bianchi et al. 2005). This is partially due to the 
fact that changing the density is equivalent to change the 
distance or luminosity, while leaving the same ionization 
parameter. 

Finally, it is important to note that, even using slightly 
different assumptions, our results showed in Fig. 1 to 
Fig. 4 are quantitatively consistent overall with those 
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showed in Fig. 1 to Fig. 6 of Bianchi et al. (2005) and 
Fig. 4 of Risaliti et al. (2005). The consistency of these 
independent studies means that the results do not dras- 
tically depend on the different assumptions and photo- 
ionization codes employed and this puts our final conclu- 
sions on a solid theoretical ground. 

3. DETAILED PHOTO-IONIZATION MODELING 
OF THE FE K ABSORPTION LINES 

In this section we describe the direct fit of the 
blueshifted Fe K absorption features detected by us in 
TA at E>6.4 keV in the spectra of the XMM-Newton 
EPIC pn with the detailed photo-ionization code Xstar. 
This allows us to derive estimates of the fundamental 
parameters of the absorbing material, such as the ion- 
ization level, column density and outflow velocity, in a 
physically self-consistent way. 

3.1. The Spectral Energy Distribution 

A good knowledge of the Spectral Energy Distribu- 
tion (SED) of each source is important for the successive 
photo-ionization modeling of the Fe K absorbers and is 
one of the main input quantities required by the Xstar 
code. The shape of the SED in different energy bands de- 
pends on the emission processes which dominate the pro- 
duction of that radiation. As stated by Elvis et al. (1994), 
a single, nearly horizontal power-law is indeed a good 
phenomenological representation of the SED of typical 
AGNs, from the IR up to the X-rays. However, some 
deviations can be observed. 

For instance, the millimeter break, which is a drop in 
the power output in the sub-millimeter band compared 
to the optical. If the drop is ~2 orders of magnitude 
the sources are classified as radio-loud, instead if higher, 
even ~5-6 decades, they are labeled as radio-quiet. The 
latter are far more common, about ~5-10 times more 
than the former (e.g., Kellerman et al. 1989 and refer- 
ences therein). In particular, the radio- loudness parame- 
ter defined as Rl = \og(fsGHz/ /4400) > 1 is usually used 
to formally differentiate between radio-loud and radio- 
quiet AGNs, where f^,GHz is the core 5 GHz radio flux 
and /4400 is the flux at 4400 A, both in units of mJy (e.g., 
Kellerman et al. 1989 and references therein). This low 
energy emission is usually associated with the power-law 
synchrotron radiation from the relativistic radio jet. 

The UV bump (or big blue bump) is a rise of the 
optical-UV continuum above the IR. It is often inter- 
preted in terms of thermal emission from the putative 
accretion disk. In particular, the beginning of the bump 
is marked by an inflection between 1 and 1.5 /im, which 
is called the near-IR inflection point. This is the only 
continuum feature which wavelength is well defined and 
it will be fundamental for the successive calculation of 
the average SED (e.g., Elvis et al. 1994). 

However, these "low energy" parts of the SED are al- 
most negligible for our purposes of the photo-ionization 
modeling of highly ionized iron ions. Instead, it is more 
important to have a good knowledge of the SED in the 
X-ray band, especially at E>1 keV. The X-ray emission 
of radio-quiet AGNs is generally supposed to come from 
inverse Compton of optical-UV seed disk photons in a hot 
corona. This component is well represented by a power- 
law continuum with typical photon index T~2, where 



FeckE r ph s 1 cm 2 keV 1 , and exponential cut-off 
at E>100 keV (e.g., Haardt & Maraschi 1991; Dadina 
2008 and references therein). 

As already stated, we need an estimate of the SED to 
calculate the Xstar grid for the successive fitting of the 
Fe K absorption lines. However, given the large number 
of sources (42) and observations (101) of the TA sample, 
it would be too much time consuming to calculate one for 
each case. Moreover, the observed flux values at different 
energies have not been measured simultaneously and the 
SEDs should be corrected for possible intrinsic source 
variability. Thus, there is the need to find a reasonable 
approximation of the typical phenomenological SED. 

We focus on the SEDs of the radio-quiet type 1 sources 
of our TA sample. This is due to the fact that possible ob- 
scuration and stellar contamination effects can prevent us 
from deriving a detailed knowledge of the intrinsic shape 
of the SED for the type 2s. We will make the reasonable 
assumption that their central engine is analogous to that 
in type Is and that the overall shape of their intrinsic 
SED is essentially the same. This is consistent with the 
widely accepted unification model for AGNs which states 
that the main difference between the two classes is essen- 
tially due to a different observing angle with respect to 
the line of sight (e.g., Antonucci 1993; Urry & Padovani 
1995 and references therein). The type Is also constitute 
the majority of the objects of our sample, 35 over 42, 
and therefore they allow for better statistics. 

We derive the average SED of the radio-quiet type 1 
sources as follows: 1) we downloaded the observed SEDs 
of the 35 type 1 AGNs from the NED databaseQ. These 
are composed by a collection of flux values measured from 
radio up to 7-rays. 2) we followed the approach of Elvis 
et al. (1994) and normalized all the SEDs to the same ar- 
bitrary flux value at the near-IR inflection point at ~1.25 
(Um (E~l eV, see above). Interestingly, the overall shape 
of the SEDs seems to be rather simple and is indeed sim- 
ilar for all the sources, see Fig. 5 (black lines). Then, 3), 
we calculated the mean flux value for each energy and de- 
rived an average SED for each object. 4) we averaged the 
flux values for each energy among all SEDs and derived 
a single average SED for all the sources. Then, 5), we in- 
terpolated the fluxes between consecutive energy points 
with power-laws. In doing so, we noted that the slopes 
were consistent one with each other for several wide en- 
ergy intervals and this suggested to us the identification 
of three main energy bands. Finally, 6), we derived that 
essentially only four energy-flux points are required to 
properly describe the overall phenomenological SED and 
we calculated the power-law slopes among them. The 
values are reported in Table 2. 

The final average phenomenological SED can be simply 
divided in three main intervals: from radio to millimetric 
(r~2), from millimetric to IR (r~0.1) and from IR to 
X-rays (r~2). We remind the reader that the relation 
between the photon index (r) and spectral index (a) is 
a = — r + 1. As shown in Fig. 5 (red line), this indeed 
provides a good representation of the overall phenomeno- 
logical SED of the type 1 AGNs. It is important to note 
that the value of T~2 for the X-ray continuum, third in- 
terval in Table 2, is in agreement with the more detailed 
work of Elvis et al. (1994) and with extensive broad-band 

2 http://nedwww.ipac.caltech.edu 
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X-ray spectral analyses (e.g., Dadina 2008 and references 
therein). Moreover, the mean phenomenological power- 
law slope in the E=4-10 keV band of the type 1 sources 
of TA is <r>~1.8. This is also expected from the two 
phase Comptonization models, which describe the pro- 
duction mechanism of the X-ray radiation (e.g., Haardt 
& Maraschi 1991), and T^2 is also usually assumed in 
photo- ionization studies of Fe K absorbers in AGNs (e.g., 
Bianchi et al. 2005; Cappi et al. 2009; TB and references 
therein) . 

3.2. Fit of the XMM-Newton EPIC pn spectra 

We consider our previous paper TA as the starting 
point for the modeling of the Fe K absorbers and we 
refer the reader to that for a detailed description of the 
data reduction, the continuum and emission lines model- 
ing and the absorption lines detection significance. Here 
we consider again the E=3.5-10.5 keV energy band of the 
28 XMM-Newton EPIC pn observations in which narrow 
Fe K absorption features at E>6.4 keV were detected 
by TA, see their Table A. 2. The spectral analysis was 
carried out using the XMM-Newton SAS 11.0 and hea- 
soft v. 6.10 packages and XSPEC v. 12.6.0. We binned 
the spectra to 1/5 of the FWHM of the EPIC pn using 
the SAS task specgroup in order to sample the resolu- 
tion of the detector. Given the X-ray brightness of the 
sources in the sample, this allows us to retain sufficient 
spectral information without the risks of oversampling 
the instrumental resolution. Given that FWHM~160- 
180 eV in the Fe K band (Guainazzi et al. 2010) this 
is equivalent to an energy binning of ~30-40 eV. How- 
ever we checked that repeating the analysis considering 
a grouping of only 25 counts per bin and the use of the 
X 2 statistic, or alternatively not grouping the data at all 
and using the C-statistic (Cash 1979), provide equiva- 
lent results. The detailed description of the fits for the 
observations of each source is reported in Appendix A. 
The errors are reported at the la level, if not otherwise 
stated. 

We started by fitting the data using the baseline mod- 
els listed in Table A. 2 of TA with the absorption lines 
modeled as inverted Gaussians and derived an estimate 
of their velocity width, o~ v . This is important for the suc- 
cessive calculation of different Xstar grids. As already 
stated in TA, the majority of the lines are unresolved 
and there we assumed widths of 10-100 eV. Given the 
limited spectral resolution of the XMM-Newton EPIC 
pn in the considered Fe K band of FWHM~160-180 eV 
(Guainazzi et al. 2010) and the S/N of the observations 
at E>7 keV, absorption lines with velocity widths lower 
than ^2000- 4000 km/s can not be resolved. In fact, re- 
peating the analysis, we could place only upper limits 
in the range <r„<3, 000-10, 000 km/s for the majority of 
them. However, in four cases we have been able to actu- 
ally resolve the width of the features to ^5,000 km/s, 
namely NGC 4151 (0402660201) cr„=5100t^ km/s, 
NGC 4051 (0157560101) <7„=4700t^$j km/s, Mrk 766 
(0304030301) cr„=4600ti^ km/s and PG 1211+143 

(0112610101) <7„=4500t^ km/s. These values are re- 
ported in Appendix A for each observation. 

The identification of spectral features in the Fe K band 
is relatively secure since only K-shell lines of heavy ions 
are expected at E>6 keV and those of iron are by far the 



most intense (e.g., Kallman et al. 2004). In TA we as- 
sumed a phenomenological identification with Fe XXV 
Hea or Fe XXVI Lya, favoring the latter if no con- 
trary indications were present, and calculated the relative 
blueshifted velocity. There we explained the reasons sup- 
porting this choice, but postponed to this second paper 
for a conclusive test. 

The possible contribution from lower ionization species 
of iron and the presence of ionized Fe K edges should be 
taken into account. These edges have rest-frame energies 
in the range from E~7.1 keV to E~9.3 keV, depending on 
the ionization state of iron, from neutral to H-like (e.g., 
Kallman et al. 2004) . As a sanity check, in TA we tested 
that the alternative modeling of the absorption lines with 
simple sharp absorption edges (zedge in XSPEC) did not 
significantly improve the spectral fits, as expected from 
the narrowness of the observed spectral features. More- 
over, it is important to note that the commonly held view 
of sharp Fe K edges is an oversimplification of the real 
process and could lead to misleading results. In fact, it 
has been demonstrated that if the adequate treatment 
of the decay pathways of resonances converging to the 
K threshold is properly taken into account, the resulting 
edges are not sharp but smeared and broadened (Palmeri 
et al. 2002; Kallman et al. 2004). This effect can be neg- 
ligible for neutral or extremely ionized iron (He/H-like) 
but is quite relevant for intermediate states, with en- 
ergies in the range E~7.2-9 keV. Furthermore, intense 
Fe K resonance absorption lines from different ioniza- 
tion states would be expected to accompany the edges. 
A proper physically self-consistent modeling of the ab- 
sorbers in the Fe K band can only be performed using a 
detailed photo-ionization code. Here we use Xstar. 

For this reason, we derived three Xstar grids consider- 
ing column densities in the range Nh =10 20 -10 24 cm~ 2 , 
ionization parameters in the interval log£=0-8 erg s _1 cm 
and turbulent velocities of 1000 km/s, 3000 km/s and 
5000 km/s, respectively. We approximated the nuclear 
X-ray ionizing continuum with the average SED de- 
rived in §3.1. In the input energy range of Xstar, from 
E=0.1 eV to E=10 6 eV, this is essentially a power-law 
with T=2 and cut-off at E>100 keV, see Table 2. We as- 
sumed standard Solar abundances (Asplund et al. 2009). 
The free parameters for each Xstar grid are the col- 
umn density Nh, logarithm of the ionization parameter 
log£ and observed absorber redshift z D . The observed 
absorber redshift is simply related to the intrinsic ab- 
sorber redshift in the source rest frame z a as (1 + z Q ) — 
(1 + z a )(l + z c ), where z c is the cosmological redshift of 
the source. The velocity can then be determined using 
the relativistic Doppler formula l + z a = (1 — /3/1 + /3) 1 / 2 , 
as (3 — v/c is here required to be positive for an outflow. 

Therefore, we fitted the data with the baseline mod- 
els listed in Table A. 2 of TA, this time without the in- 
clusion of the inverted Gaussian lines, and modeled the 
absorption features using Xstar. For each observation, 
if the velocity width cr v of the lines was constrained, we 
performed a fit using the Xstar grid with the closest tur- 
bulent velocity, which turned out to be possible only for 
the four cases reported above with o~ v ~5,000 km/s. In- 
stead, if only an upper limit was derived, we performed 
fits using different Xstar grids with turbulent velocities 
of 1000 km/s, 3000 km/s and 5000 km/s until these val- 
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ues were lower or equal to the line velocity width and 
then checked which one provides a better fit. 

To find the Xstar solution(s), i.e. the possible minima 
in the x 2 distribution, we stepped through the absorber 
rcdshift in small increments of Az=10~ 3 in the interval 
between 0.1 to —0.4, leaving also the other parameters 
of the absorber, i.e. Nh and log£, and the continuum 
free to vary and recorded the x 2 values. For the case 
of 1H0419-577 we stepped the redshift starting from 0.15 
because the source cosmological redshift is £=0.104. This 
allows us to blindly sample the whole 6-10 keV energy 
band and search for possible alternative Xstar solutions. 
The x 2 °f the fits plotted against the absorber redshift 
is reported for all observations in Fig. 12 in Appendix B. 
This is the best way to derive a physically self-consistent 
identification of the absorption features. In fact, Xstar 
simultaneously takes into account all the resonance ab- 
sorption lines and edges from all elements up to Z=30 
for different ionization states. The detailed description 
of the Xstar fits for each observation is discussed in Ap- 
pendix A and the summary of the best-fit Xstar param- 
eters is reported in Table 3. 

The possible presence of different minima in the x 2 dis- 
tribution indicates that more Xstar solutions are possible 
for distinct rcdshifts or turbulent velocities. In particu- 
lar, there are 6/28 cases with two redshift solutions with 
equivalent x 2 levels, indicating a possible degeneracy in 
the identification of the absorption line(s) as resonances 
from Fe XXV or Fc XXVI. In these cases, we consid- 
ered the best-fit model to be the one associated with the 
lowest x 2 value and the relative fit statistic and F-test 
probability are reported in Table 3. Then, if the possible 
alternative solutions had x 2 values consistent at the 90% 
level with respect to the best-fit, i.e. with a difference 
of Ax 2 <2.7, we averaged the values of the parameters 
among them and derived the errors as half of the inter- 
val between the lowest and the highest possible values. 
This allows us to take into account both the statistical 
errors of measure and the uncertainty in the determi- 
nation of the best-fit model or, equivalently, of the line 
identification. However, it is important to note that all 
the detected narrow Fe K lines are indeed consistent with 
an identification as resonance transitions essentially from 
Fe XXV and Fe XXVI ions. This is also demonstrated 
by the overall consistency of the outflow velocities of the 
absorbers reported in Table 3 and in Table A. 2 of TA. 
We note that in all the observations only one single Xs- 
tar component was required to properly model the Fe K 
absorption lines. 

In Fig. 6, Fig. 7 and Fig. 8 we report the source and 
background spectra and the best-fit model for three rep- 
resentative cases. More details on the spectral fitting 
and the best-fit parameters are reported in Appendix A 
and Table 3 for all the observations, respectively. Fig. 6 
shows the case of PG 1211+143 (0112610101) in which 
a UFO with iw=45,300 ± 900 km/s (0.151 ± 0.003c) 
has been detected with a high confidence of Pf>99.99%. 
The detection in this observation was already reported 
by Pounds et al. (2003). The line is resolved and the ve- 
locity width is cr„=4, 500^200 km/s. Thus, the absorber 
best-fit is provided by an Xstar grid with turbulent ve- 
locity of 5,000 km/s. As it can be seen from Fig. 12 in 
Appendix B, there is only one solution, which indicates 



the absorption is mainly due to Fe XXV and intermediate 
states Fe XX-XXIV transitions. The estimated column 
density and ionization are -/V#=(8.0t 2 ' 2 ) x 10 22 cm~ 2 

and log£=2. 87^0^0 er g s_1 cn b respectively. These val- 
ues are consistent with the results in Fig. 2. In fact, 
we can note that the high associated absorption line 
EW of -130 eV (see Table A.2 of TA) for the Fc XXV 
Hea transition can be reached with an ionization of 
log£~2.9 erg s -1 cm for iV#~10 23 cm~ 2 and a velocity 
broadening of ^5,000 km/s. From Fig. 3 we see that for 
the estimated ionization parameter and velocity broad- 
ening, the Fe XXVI Lya is ~10 times weaker than the 
Fe XXV Hea and this explains why it is not observable 
in the spectrum. Finally, from Fig. 4 we can see that for 
a velocity width of ^5,000 km/s, the Fe K lines are not 
saturated up to column densities of ^10 23 -10 cm -2 . 
Due to the relatively low ionization level, in the best-fit 
model in Fig. 6 we can note the presence of additional 
weak absorption structures at high energies due to the 
several resonance series converging to the Fe XX-XXV 
edges (e.g., Kallman et al. 2004). 

Fig. 7 shows the case of Mrk 79 (0400070201) in which 
we detected a UFO with v out =27, 600 ± 1,200 km/s 
(0.092 ± 0.004c) with a high confidence of P F =99.97%. 
The relative absorption line is not resolved, but we can 
place an upper limit of a v < 6,400 km/s. From Fig. 12 
we can note that there is only one redshift solution, 
consistent with an identification as mainly absorption 
from Fc XXVI. As discussed in Appendix A, the best 
fit is provided by an Xstar grid with turbulent veloc- 
ity of 1,000 km/s (see middle panel in Fig. 7) but the 
one with 3,000 km/s is equivalent at 90%. The column 
density is 7Vh=(19. 4 ± 12.0) x 10 22 cm~ 2 and ioniza- 
tion log£=4.19 ± 0.23 erg s" 1 cm. The EW of the as- 
sociated blue-shifted absorption line detected by TA is 
^40cV. Consistently with Fig. 2, this can be reached 
with 7V#~10 23 cm~ 2 and log£~4 erg s _1 cm for a ve- 
locity width of ^1,000-3,000 km/s. From the ratios in 
Fig. 3 we note that for these values the Fc XXVI Lya 
is indeed more intense, ~5-10 times, than the Fe XXV 
Hea and this is consistent with the fact that only one 
absorption line could be detected at high significance 
given the limited S/N of the spectrum. Moreover, from 
the ratios in Fig. 4 we see that for velocity widths of 
~1, 000-3,000 km/s, the ls-2p lines start to be saturated 
at AT#~10 23 cm -2 , instead the other ls-3p transitions 
at higher energies can still grow and their intensities are 
a bit higher than expected from the simple ratio between 
the relative oscillator strengths. 

Finally, in Fig. 8 we show the case of NGC 3516 
(0401210601) in which we detected a series of lines in- 
dicating the presence of a lower velocity system with 
v out =3, 300 ± 600 km/s (0.011 ± 0.002c) with a high con- 
fidence of Pp>99.99%. The complex absorption in this 
source was already discussed by Turner et al. (2008). 
The lines are unresolved, but the upper limit on the ve- 
locity width is a t ,<2,600 km/s. From Fig. 12 and Ap- 
pendix A we can note only one best-fit solution consistent 
with absorption mainly from Fe XXV-XXVI. The best- 
fit is provided by an Xstar grid with turbulent velocity 
fo 1,000 km/s. The column density is N H = (6.7t 2 ; 3 ) x 
10 22 cm~ 2 and ionization log£=4.05l ; 08 erg s" 1 cm. 
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The EW of the associated four absorption lines detected 
by TA is ~25-30 eV and the estimated parameters are 
quantitatively consistent with what is shown in Fig. 2. 
Given these parameters, we can see from Fig. 3 and Fig. 4 
that the ls-2p lines are slightly saturated and their ra- 
tios with respect to the ls-3p is slightly lower than what 
simply expected from their relative oscillator strengths. 

3.3. Consistency checks of the spectral analysis and 
modeling 

The detailed treatment of the continuum modeling 
complexities, the background checks, the consistency 
check with the simultaneous XMM-Newton EPIC Metal 
Oxide Semi-conductor (MOS) camera and the detection 
significance of the lines through the F-test and exten- 
sive Monte Carlo simulations were already reported in 
TA. As shown in Fig. 6, Fig. 7 and Fig. 8, here we note 
again that the background level of the observations of 
the sample in the 7-10 keV band is always less than 
~10% of the source counts. We also checked that none 
of the reported blue-shifted absorption lines were induced 
by an erroneous subtraction of the EPIC pn instrumen- 
tal background emission lines from Ni Ken and Cu Ka 
at the energies of 7.48 keV and 8.05 keV, respectively 
(Katayama et al. 2004). In fact, only 4/22 of the ab- 
sorption lines detected at E>7 keV may have observed 
energies consistent at the 90% level with these and the 
features are present in the source spectra with and with- 
out background subtraction. Moreover, the detected fea- 
tures can not be attributed to some sort of EPIC pn cal- 
ibration artifact because they have been detected always 
at different energies and are narrow. Instead, if due to 
instrumental calibration problems, they would have been 
expected to be present always at some specific energies 
or, if due to effective area calibration problems, to induce 
broad continuum distortions. 

Here, we perform additional tests to check the con- 
sistency of the photo-ionization modeling. We start with 
the SED. We checked that a different choice of the power- 
law slope in the wide range r=1.5-2.5 has <10% effects 
on the Nh and log£ estimates, within the measurement 
errors. This is consistent with what derived in §2.2. We 
then compared the results using the average phenomcno- 
logical SED derived in §3.1 with the power-law of T=2 
assumed in the curve of growth analysis of §2.2. We 
reproduced the EW distributions reported in Fig. 2 for 
these two cases and find a slight shift toward higher ion- 
ization parameters for the average SED to obtain the 
same value of the power-law SED at the peak. This ef- 
fect is essentially due to the different high energy part 
of these two SEDs. In fact, the power-law and average 
SEDs have the same slope of T=2 until E=10 5 eV, then 
the average SED drops to zero and the power-law one 
instead extends up to the upper limit of the Xstar in- 
put energy band of E=10 6 eV. This means that, given 
the same bolometric luminosity, the average SED has 
slightly less high energy photons that can effectively ion- 
ize iron with respect to the power-law one. Therefore, to 
obtain the same iron ionic populations, the average SED 
needs to increase the global ionization level of the ma- 
terial. However, the presence of the high energy cut-off 
at E~100 keV in the average SED causes only a neg- 
ligible difference in the estimated ionization parameters 
in Table 3, being less than 5% and within the measure- 



ment errors. Moreover, this high energy cut-off is natu- 
rally expected from comptonization models of the X-ray 
continuum in AGNs (e.g., Haardt & Maraschi 1991 and 
references therein) and is in agreement with what is ob- 
served for local Seyfert galaxies (e.g., Dadina 2008 and 
references therein). Moreover, as already noted by McK- 
ernan et al. (2003), the presence or absence of the possi- 
ble UV-bump in the SED has a negligible effect on the 
parameters of the photo-ionized gas in the Fe K band. 
Again, this and the previous results are essentially due to 
the fact that only the high energy part of the X-ray con- 
tinuum at E>6 keV is capable of photoionizing the highly 
ionized absorbers studied here, which are dominated by 
Fe XXV/XXVI. Therefore, we can directly compare the 
Xstar results in Table 3 with the general trends reported 
in the curve of growth analysis of §2. 

Another important parameter of the Xstar grids is the 
turbulent velocity of the plasma. It should be noted that 
the turbulent velocity considered in the grid is different 
from the velocity broadening parameter b considered in 
the previous curve of growth analysis (§2.2), which is the 
square root of the quadratic sum of the turbulent and 
thermal velocities. Xstar already incorporates the ther- 
mal broadening of the lines. However, as already stated 
in §2.2, the thermal velocity term for the highly ionized 
absorbers considered here is essentially negligible and the 
velocity broadening parameter b is dominated by the tur- 
bulent component. Therefore, for the highly ionized Fe K 
features studied here and if not strongly affected by line 
blending, the velocity parameter b defined in §2.2, the 
turbulent velocity of the Xstar grids Vt and the velocity 
width of the lines a v measured from the spectra are all 
comparable. We tested that for lower/higher choices of 
this parameter, in the wide range ^100-5000 km/s, the 
resultant estimate of the ionization parameter is essen- 
tially not affected, although the derived absorber column 
density is found to be slightly higher/lower. This is due 
to the line saturation effect, which changes the slope of 
the curve of growth for high column densities, see Fig. 1. 
To compensate for these possible systematics, we used an 
Xstar grid with a turbulent velocity closest to the line 
velocity width, for the few cases in which the latter was 
constrained from the data. Instead, for the other cases 
with only upper limits, we averaged the values of the 
parameters derived from fits using grids with different 
turbulent velocities, see §3.2 and Appendix A. 

4. RESULTS 

The summary of the best-fit parameters of our photo- 
ionization modeling of the Fe K absorbers detected in 
the XMM-Newton EPIC pn spectra of TA is reported in 
Table 3. We can see that the absorbers are systemati- 
cally outflowing and their velocities are determined with 
high accuracy, with errors taking into account also the 
uncertainty on the line identification (sec § 3.2 and Ap- 
pendix A). The velocity values are consistent overall with 
those derived in TA using inverted Gaussians to model 
the absorption lines and a more simplistic line identifica- 
tion, see §4.1 and Table A. 2 of TA. The observations in 
which we detected a UFO are marked with an asterisk 
in Table 3. 

The fraction of sources showing UFOs with respect to 
the whole sample is again 15/42, which corresponds to 
^35%. However, as already discussed in §4.2 of TA, if 
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we take into account the fact that we have only con- 
servative estimates on the velocities and that there is a 
number of spectra in the sample with low S/N, the frac- 
tion of UFOs is more probably in the range — 40%-60%. 
The fraction of sources in the sample with UFOs having 
velocities >0.1c is —11/42 (—25%). This frequency is in- 
stead —11/15 (—70%) only considering those with UFOs. 
The fraction of sources showing UFOs with respect to 
the total number with detected Fe K absorbers is 15/19 
(—80%). Therefore, the great majority of sources with 
detected Fe K absorbers show the presence of UFOs. The 
number of observations with detected Fe K absorption 
lines consistent with UFOs is 19/28 (—70%), as reported 
in Table 3. 

The histograms representing the distribution of the 
outflow velocity, ionization and column density for the 
UFOs and the lower velocity systems are shown in Fig. 9, 
Fig. 10 and Fig. 11, respectively. Following the same 
considerations of TA, we averaged the values reported in 
Table 3 among the observations of each source and con- 
sidered the fractions with respect to the total number of 
sources in the sample (42). Instead, using each observa- 
tion separately would have introduced a bias because the 
sources with more observations would have had a higher 
weight. From Fig. 9 we can see that the velocity dis- 
tribution of UFOs spans from -10,000 km/s (-o!033c) 
up to '-TOO, 000 km/s (—0.3c), and has a peak and mean 
value at —42,000 km/s (— 0.14c). Instead, by definition, 
the distribution of the lower velocity systems is narrow 
and has a mean value of —3,000 km/s (—0.01c). 

However, it should be noted that the present data are 
affected by several selection effects that hampers the de- 
tection of absorption lines at higher energies, or equiv- 
alently, of absorbers with higher outflow velocities. An 
estimate of the magnitude of this effect was already re- 
ported in §4.2 and Fig. 6 of TA. This is due to a combi- 
nation of different causes, such as: the upper limit of the 
EPIC pn detector bandwidth at —10 keV, which does not 
allow the detection of Fe K lines with blueshifts higher 
than —0.3-0. 4c for the local sources studied here, the 
drop in effective area and energy resolution of the instru- 
ment at E>7 keV, the intrinsic lower number of counts 
for the power-law spectra at high energies and the rel- 
atively short exposures. We can derive an estimate of 
the influence of the limited S/N of the XMM- Newton 
observations on the derived UFO velocity distribution 
normalizing the fraction of observations for each velocity 
interval in the histogram taking into account the number 
of observations that actually offer enough statistics to de- 
tect an absorption line at > 3a. As already discussed in 
§4.2 of TA, this number decreases rapidly for lines with 
lower EW and higher associated velocity. Consequently, 
the allowed fraction of observations with higher velocity 
increases. Assuming the typical absorption line EWs of 
20-50 eV, this implies that the peak and mean of the 
outflow velocity distribution of UFOs in Fig. 9 can po- 
tentially shift from —0.14c up to — 0.16-0. 18c. 

As reported in Table 3, we have been able to mea- 
sure the ionization parameter of the absorbers with good 
accuracy for all observations with Fe K absorbers, tak- 
ing also into account the uncertainty in the line iden- 
tification. The distributions of the ionization param- 
eter of the Fe K absorbers ascribable to UFO and to 
the lower velocity systems are shown in Fig. 10. The 



values log£— 2.5 and 6 erg s 1 cm delimit the range in 
which Fe XXV '/XXVI lines have observable EW>10 eV, 
consistently with Fig. 2. We can note that the ioniza- 
tion distribution for the UFOs is much broader than for 
the other lower velocity absorbers, spanning the whole 
range from log£ of 2.5 up to 6 erg s _1 cm, with a mean 
value of log£— 4.2 erg s _1 cm. The very high log£ values 
of the UFOs may possibly suggest that in some cases 
the material could not be detected in absorption due 
to the fact that even iron was completely ionized. In- 
stead, the lower velocity Fe K absorbers are limited in 
the range log£— 2.5-4.5 erg s _1 cm, with a mean value of 
log£— 3.5 erg s" 1 cm. 

Fig. 1 1 shows the distribution of the column densities 
for the UFOs and the lower velocity systems. We in- 
cluded only the best-fit values in Table 3 for the cases 
in which the Nh was constrained. We limited the treat- 
ment to Compton-thin absorbers with iV#<10 24 cm~ 2 . 
We can see that both distributions have column densi- 
ties in the wide range iVf/— 10 22 -10 24 cm" 2 , with those 
of the UFOs slightly higher than the lower velocity sys- 
tems. Their mean values are of A?#— 1.3 x 10 23 cm -2 and 
iVf/— 6.5 x 10 22 cm~ 2 , respectively. This might suggest 
the presence of possible Compton thick material, at least 
in the inner regions of the flow. 

Finally, it should be noted that the velocities and col- 
umn densities of the Fe K absorbers derived by fitting the 
X-ray spectra depend on the unknown inclination angle 
of the outflows with respect to our line of sight. In other 
words, they depend whether we are actually looking di- 
rectly down the outflowing stream or if we intercept only 
part of it (e.g., Elvis 2000). Therefore, these measures 
should be regarded only as conservative estimates and 
the intrinsic values might be higher. 

5. CONCLUSIONS 

In this paper we performed an extension of the work 
presented by us in TA, in which we showed the detection 
of blueshifted Fe K absorption lines ascribable to UFOs 
in >35% of 42 local radio-quiet AGNs observed with 
XMM-Newton. Here, we started with a curve of growth 
analysis of the main Fe XXV/XXVI K-shell transitions 
in photo-ionized plasmas. We studied the dependence 
of the equivalent width of the absorption lines as a func- 
tion of the ionization parameter, column density, velocity 
broadening and the lines ratios. Then, we estimated an 
average SED for the sample sources and directly mod- 
eled the Fc K absorbers in the XMM-Newton EPIC pn 
spectra with the photo-ionization code Xstar. A single 
absorption component is capable of self consistently pa- 
rameterizing the simple Gaussian absorption lines of TA 
for each observation. We derived the absorbers column 
density, ionization parameter and outflow velocity. The 
summary of the best-fit parameters is reported in Ta- 
ble 3. 

We confirm that the frequency of sources showing 
UFOs in the radio-quiet sample is >35% and that the 
majority of the Fe K absorbers are indeed associated 
with UFOs (>80%). We derived that the outflow veloc- 
ity distribution of the UFOs spans from — 10,000 km/s 
(~0.03c) up to -100,000 km/s (-0.3c), with a peak and 
mean value at -42,000-55,000 km/s (-0.14-0.18c), tak- 
ing into account also the presence of selection effects that 
hampers the detection of the higher velocity absorbers. 
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These effects can only be reduced by performing much 
longer observations of these sources with the present in- 
struments or using the higher sensitivity and energy res- 
olution of the calorimeters on board the future Astro-H 
and the proposed Athena missions. 

The ionization parameter of the UFOs is very high and 
in the wide range log£~3-6 erg s _1 cm, consistent with 
Fe XXV/XXVI being the dominant lines, with a mean 
value of log£~4.2 erg s _1 cm. This might give rise to the 
possibility of the presence of additional columns of ma- 
terial, which are unseen to absorption line spectroscopy 
because even iron is completely ionized. The column 
densities associated with the UFOs are also large and in 
the broad range 7Vy~10 22 -10 24 cm~ 2 , with a mean value 
of A/#~10 23 cm -2 . This might suggest the presence of 
Compton thick material, at least in the inner regions of 
the flow. However, besides possible selection effects, the 
measured velocities and column densities should also be 
regarded as conservative estimates because they depend 
on the unknown inclination angle of the flow with respect 
to the line of sight. 

Depending on the actual outflow geometry and cover- 
ing factor, the latter estimated from the detection fre- 
quency on the whole sample of ~0.4-0.6 in §5 of TA, 
this material might be observed also in emission, through 
fluorescence/recombination lines, continuum reflection, 
bremsstrahlung and/or black body emission. The elec- 
tron temperature of kT~100 eV (see §2.2) then suggests 
the continuum emission might be mainly in the soft X- 
rays (e.g., King & Pounds 2003; Pounds & Reeves 2007; 
Sim et al. 2008; Sim et al. 2010 and references therein). 
Finally, it is important to note that the parameters here 
derived for the UFOs in the local radio-quiet AGN sam- 
ple of TA in Table 3 are consistent overall with those 
estimated for the complementary, but much smaller, lo- 
cal radio-loud AGN sample in TB, suggesting a possible 
link between the two AGN classes. 

We refer to a successive paper III of this series 
(Tombesi et al. in prep.) for a detailed discussion of the 
physical properties, geometry and energetics of the UFOs 
here analyzed and also a comparison with the classical 
AGN warm absorbers. Finally, significant improvements 
are expected from the high effective area and supreme 
energy resolution in the Fe K band offered by the X-ray 
microcalorimeters on board the future Astro-H and the 
proposed Athena missions. 
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1Q20 10 21 10 22 10 23 10 24 



N H (cm- 2 ) 

Figure 1. Curve of growth for the Fe XXV Hect (r + i) (upper panel) and Fe XXVI Lya (ri + T2) (lower panel) lines, for different values 
of the velocity broadening parameter b. From bottom to top b=20, 100, 500, 1000, 3000, 5000, 10000 km/s. The ionization parameter is 
fixed to log£=3.5 erg s _1 cm. 
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Figure 2. Estimated EW of the Fe XXV Hea (r + i) (upper panel) and Fe XXVI Lya (ri + rz) (lower panel) transition as a function of 
log£, for different column densities. The velocity broadening is 6=100 (black), 1000 (red), 5000 km/s (blue). 




Figure 3. Ratio of the EWs of the Fe XXV Hea (r + i) and Fe XXVI Lya (ri + 1:2) lines as a function of log£, for different column 
densities. The velocity broadening is b= 100 (black), 1000 (red), 5000 km/s (blue). 
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Figure 4. Ratio of the EWs of the Fe XXV Hea/He/3 (r + i) (upper panel) and Fe XXVI Lyoj/Ly/3 (ri + 12) (lower panel) lines as a 
function of the column density, for different values of the velocity parameter b. The ionization parameter is fixed to log£=3.5 erg s" 1 cm. 
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Figure 6. XMM-Newton EPIC pn spectrum of PG 1211+143 (0112610101) in the 4-10 keV band in which a UFO with iw~0.15c has 
been detected. Upper panel: background subtracted source spectrum with superimposed the best-fit model. Middle panel: best-fit model 
using an Xstar grid with turbulent velocity of 5000 km/s. Lower panel: background spectrum. Data are binned to 1/5 of the FWHM of 
the instrument. 
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Figure 7. XMM-Newton EPIC pn spectrum of Mrk 79 (0400070201) in the 4-10 keV band in which a UFO with iw~0.1c has been 
detected. Upper panel: background subtracted source spectrum with superimposed the best-fit model. Middle panel: best-fit model using 
an Xstar grid with turbulent velocity of 1000 km/s. Lower panel: background spectrum. Data are binned to 1/5 of the FWHM of the 
instrument. 
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Figure 8. XMM-Newton EPIC pn spectrum of NGC 3516 (0401210601) in the 4-10 keV band in which a lower velocity system with 
i> ou t~0.01c has been detected. Upper panel: background subtracted source spectrum with superimposed the best-fit model. Middle panel: 
best-fit model using an Xstar grid with turbulent velocity of 1000 km/s. Lower panel: background spectrum. Data are binned to 1/5 of 
the FWHM of the instrument. 
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Figure 9. Histograms representing the distribution of the outflow velocity of UFOs (solid line) and the lower velocity systems (dashed 
line), normalized to the total number of sources in the sample. 
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Figure 10. Histograms representing the distribution of the ionization parameter of the as UFOs (solid line) and the lower velocity systems 
(dashed line), normalized to the total number of sources in the sample. 
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Figure 11. Histograms representing the distribution of the column density of the UFOs (solid line) and lower velocity systems (dashed 
line), normalized to the total number of sources in the sample. 
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Table 1 

Parameters of the considered Fe XXV and Fc XXVI K-shcll transitions. 



Ion 


ID 


Transition 


<E> 


Line 


E 


flu 


A u i 








koV 




(kcV) 


(xicr 2 ) 


(xlO 12 ) 


Fc XXV 


Hea 


ls 2 -ls2p 


6.697 


(r) 
(i) 


6.700 


70.4 


457 










6.668 


6.9 


44 




He/3 


ls 2 -ls3p 


7.880 


(r) 


7.881 


13.8 


124 










(i) 


7.872 


1.7 


1 


Fc XXVI 


Lya 


ls-2p 


6.966 


(n) 


6.973 


28.0 


296 










(ra) 


6.952 


14.0 


293 




Ly/3 


ls-3p 


8.250 


(ri) 


8.253 


5.3 


79 










(ra) 


8.246 


2.6 


78 



Note. — The line parameters have been taken from the accurate NIST atomic database. (1) Iron ion. (2) Line identification. (3) Atomic 
transition. (4) Average line energy, weighted for the oscillator strengths of its components. (5) Line components. (6) Energy of the specific 
transition. (7) Oscillator strength. (8) Einstein coefficient. 
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Table 2 

Average phenomenological SED for the type 1 objects of the TA radio-quiet AGN sample. 



Interval E (eV) T a 

1 5 X 10~ 6 -5 X 1(T 5 2 -1 

2 5 x l(T 5 -2 x 1(T 2 0.1 0.9 

3 2 x 10- 2 -l x 10+ 5 2 -1 
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Table 3 

Summary of the best-fit Xstar parameters. 



Source 


OBSID 




Vout 


log? 


N„ 


A x 2 /Au 




Pf 


INLtO 41 01 


U4U200U2U1 


tr 1 +1.8 

5-l_!. 4 


n infi+°- ()()7 

°- 106 -0.007 


4 41 +0.92 
4 - 41 -0.08 


> 2 


21.4/ o 


219.D/ 10 ( 


no ono/ 

yy.ou/o 


T/"l a ion A 


m A ^ A A f\~\ c\-\ * 

0147440101 


< 3.1 


n noa+0-004 
°- 098 -0.004 


t- o,+0.94 
°-"-'*-0.94 


> 2 


12.5/3 


240.D/159 


n"7 nnO/ 

y 7. oo 7o 


l\T/^t/~1 Q70Q 


0112210101 


< 2.6 


o.oio_ 005 


9 -, ,+0.15 
z - ±o -0.09 


9 O+0.9 
z -°-0.5 


1 1 Aj 6 


om e /icn 
2U1.D/ loU 


nn i^noy 

yy.bo/o 




ni 1 ooi noni 
0112210201 


< 2.7 


< 0.007 


9 07+O.O6 
AS '-0.05 


o Q+0-2 
z -°-0.2 


ACl All 

4U.4/0 


17Q r /1 C7 

1 io.o/ 157 


>yy.yy /o 




ni 1 ooi n cm 
0112210001 


< 2.6 


< 0.007 


9 O.K+ - 06 
z - ao -0.06 


9 7+0.2 
2 "-0.2 


77.8/3 


i no o /i C7 

lyo.o/ 15 ( 


\ nn nno/ 

>yy.yy /o 


iNLrC oolo 


0401210401 


< 2.7 


n nnK+0.002 
°- 006 -0.002 


4 - 04 -0.14 


r ,,+2.0 


58.4/3 


oner a 1 1 nT 

20o.4/iy7 


~ nn nnQ/ 

>yy.997o 




n/im oi n cm 
0401210001 


< 2.8 


0.008_ 002 


,, 1 o+0.09 
4 ' 18 -0.06 


7 n+3-2 


45.3/3 


1 nr a /1 t(? 

lyo.u/ loo 


-~ nn nnoy 

>yy.yy /o 




0401210601 


< 2.6 


n m i +0.002 

o.oni 002 


a nK+0.09 
4 - 05 -0.08 


6.7±f 3 4 


54.8/3 


170.3/154 


- nn nnOy 

>99.997o 




n a ni oi 1 nm 
0401211U01 


< 2.8 


n m ,+0.003 

0.013_ 003 


4 -, O+0.12 

4 - 18 -0.10 


7 1+2.4 
7 -l-2.2 


25.1/3 


1 ;y.D/ 154 


nn nooy 

yy. yo/o 


Mrk oUy 


m onvom m * 
OloO 7 2U1U1 


< 3.8 


n 17 o+0.004 
u -! ,,:> -0.004 


t- nfi +0.58 


> 6 


11.8/3 


i on o /i en 

ioy.2/ loU 


nn a noy 

yy. 40/0 




UoUDUyU2Ul 


< o.z 


n iQQ+0.004 

U.138_ 004 


c n9 +0.98 
°- uz -0.98 


> 


i n 7 /q 
12. / /o 


1 77 fi/1 co 

1 i ( .U/ lOo 


no nn<V 

yy.uu/o 




UoUoUyU4Ul 


< 5.0 


n 91 v +0.024 
°- 217 -0.024 


A 9C - + 1.46 
4 ' 25 -1.46 


> 0.7 


14.5/3 


om n / i cn 
2U1.U/ ioy 


nn nno/ 
yy. OO/o 


Ark 120 


m /i7inni m * 
014 / lyUlUl 


< 3.3 


n ooy+0.022 
u - z °' -0.022 


, t-c + 1.29 
4.55_! 29 


> 0.7 


12.5/3 


T /in o /i co 

149.3/153 


nn a r\(y/ 

99.40/0 


TV /T*. 1 OTn 

Mrk 279 


nono a on cn 1 

0002480501 


< 4.5 


< 0.007 


, 7 Q+0.59 
— 0.59 


A fi+3-6 
4 - 6 -3.6 


14.0/3 


TOO C /1CO 

lo2.o/152 


nn o c 0/ 

yy.oo/o 


Mrk /y 


04000/0201 


^ P. A 

< 0.4 


n nQO+ 004 
0.092_ 004 


, -iq+0.23 
±a -0.23 


1Q - + 12. 
±».l _ 12 Q 


in 1 /Q 

iy. i / o 


i qi n / i oo 

loi.y/ ioo 


no nW 

yy.y ( /o 


NGC 4051 


0109141401* 


< 4.2 


r> ^0^7+0.025 
°- 037 -0.025 


4 - 37 -1.54 


> 0.8 


13.8/3 


160.5/155 


99.50% 




uio * oouiui 


A 7+1-8 
4 ' 7 -1.5 


n 9n9 +0.006 
u.zuz_ Q 00g 


3.06_ 19 


c o+l.l 
5 - 2 -1.3 


22. U/ o 


102.D/ 101 


no nci'V 

yy.yo /o 


Mrk /oo 


nin/inQnom * 
0o040o0o01 


,1 fi+2-7 
4 - 6 -1.8 


n nao+0-006 
U.0»2_ 0Q6 


3.46_ 032 


9 i +0.4 
^•-"--0.3 


22.4/3 


lo2.o/ 154 


nn n c 07 

yy.yo /o 




nin.ininr'n-i 

0304030501 


< 4.4 


n noo+0.002 

u.o»»_ 002 


„ oo+O.OS 
4 - 28 -0.11 


in n+6-4 


21.2/3 


187.1/151 


nn nnOy 

99.907o 


Mrk 841 


0205340401* 


< 4.0 






> 1 


12.7/3 


134.9/127 


99.07% 


ESO 323-G77 


0300240501 


< 3.7 


012+ ' 004 


3.74+°; 3 37 


9 o 1 + 12.3 


77.0/3 


146.8/132 


>99.99% 


1H0419-577 


0148000201* 


< 7.0 


n n 7Q+ - 007 


3-69l°;f 3 


17 o+14.3 


14.8/3 


113.9/102 


99.43% 


Mrk 290 


0400360601* 


< 4.0 


163+ 024 

U.100_g Q 24 


•j-yi-i .17 


25 9+ 24 - 


13.5/3 


90.2/101 


99.74% 


Mrk 205 


0124110101* 


< 3.1 


100+ ' 004 


4 92+0-56 
^• az -0.56 


> 14.5 


11.9/3 


123.8/118 


98.80% 


PG 1211+143 


0112610101* 


4 -°-1.2 


n 151 +0.003 


9 oy+0.12 
z -°'-0.10 


8.0+j' 2 


70.5/3 


125.5/114 


>99.99% 


MCG-5-23-16 


0302850201* 


< 13.0 


11fi+ - 004 


4 Qo+0.08 
^•°°-0.08 




41.3/3 


257.1/158 


>99.99% 


NGC 4507 


0006220201* 


< 4.2 


n iqq+0-024 

o.iyy_ 24 


4.53li:lS 


> 0.9 


10.1/3 


148.7/144 


97.70% 


NGC 7582 


0112310201* 


< 11.0 


285+ ' 003 

U.ZOJ_g QQ 3 


O QQ + 0.09 


23.41& 1 


20.9/3 


111.0/100 


99.94% 



Note. — The errors are at the la level. The lower limits on iVjj and upper limits on v out arc at the 90% confidence level. The * marks the 
observations with detected UFO. (1) Source name. (2) XMM-Newton observation ID. (3) Absorption lines velocity width in units of 10 3 km/s. (4) 
Outflow velocity in units of c. (5) Logarithm of the ionization parameter in units of erg s _1 cm. (6) Absorber column density in units of 10 22 cm -2 . 
(7) X 2 improvement after including the Xstar table in the baseline model and additional model parameters. (8) Final best-fit % 2 after including 
the Xstar table and degrees of freedom. (9) F-test probability. 
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APPENDIX 
NOTES ON SINGLE SOURCES 
NGC 4151 

Observation 0402660201: the velocity width of the absorption line is measured as a v = 5100^140° km/s. Fig. 12 
in Appendix B shows the fit solution for the Xstar grid with v t = 5000 km/s. The parameters are log£=4.41+QQ 2 s 
erg s" 1 cm, Nh>2 x 10 22 cm -2 and rcdshift z=— 0.099 ± 0.007, which corresponds to ?w=0. 106 ± 0.007c. The initial 
baseline model statistic is x 2 =241.0 for d.o. /.=160 and the fit improvement is A% 2 =21.4 for three model parameters. 
This provides the best-fit solution, with F-test probability is P F =99.8%. We detected for the first time a UFO in this 
observation. 



IG4329A 

Observation 0147440101: the velocity width of the absorption line is a v < 3100 km/s. From Fig. 12 in Appendix B 
there is one redshift solution. The initial baseline model statistic is x 2 =253.1 for do./.=162. Using an Xstar grid with 
u t =1000 km/s we obtain (1) log£=5.17+°;77 erg s" 1 cm, N H >2 x 10 22 cm" 2 and z=-0. 078 ± 0.003, which corresponds 
to v ottt =0.097± 0.003c. The fit improvement is Ax 2 =12.5 for three model parameters. Instead, using an Xstar grid 
with v t =3000 km/s we obtain (2) log£=5.63t?;i5 er S s_1 cm > n h>2 x 10 22 cm" 2 , z=-0.079±0.003, which corresponds 
to zw=0.098± 0.003c. The fit improvement is Ax 2 =10.8 for three model parameters. The best-fit is provided by (1), 
with F-test probability Pf—97%. The fits (1) and (2) are equivalent at 90%. The detection of a possible UFO in this 
observation was already reported by Markowitz et al. (2006). 

NGC 3783 

Observation 0112210101: the velocity width of the absorption line is a v < 2600 km/s. From Fig. 12 in Appendix B, 
there is only one redshift solution. The initial baseline model statistic is x 2 =219.0 for d.o. f. =163. Using an Xstar grid 
with w t =1000 km/s we obtain log^=2.13^;J| erg s" 1 cm, iVj^ =(2.3±g;|) x 10 22 cm" 2 and z=-0.001 ± 0.005, which 
corresponds to v otl t=0.010± 0.005c. The fit improvement is Ax 2 =17.4 for three model parameters. This provides the 
best-fit, with F-test probability P F =99.6%. 

Observation 0112210201: the velocity width of the absorption lines is a v < 2700 km/s. From Fig. 12 in Appendix B, 
there is only one redshift solution. The initial baseline model statistic is x 2 — 213.9 for d.o. f. =160. Using an Xstar grid 
with u t =1000 km/s we obtain log£=2.87t£$? erg s" 1 cm, N H ={2.8±%%j x 10 22 cm" 2 and redshift z=0.009 ± 0.004, 
which corresponds to u out <0.007c. Ax 2 =40.4 for three model parameters. This provides the best-fit, with F-test 
probability P F >99.99%. 

Observation 0112210501: the velocity width of the absorption line is a v < 2600 km/s. From Fig. 12 in Appendix 
B, there is only one redshift solution. The initial baseline model statistic is x 2 =271.6 for do./.=160. Using an Xstar 
grid with w t =1000 km/s we obtain log£=2.98±^6 erg s" 1 cm, N H =(2.7±%%) x 10 22 cm" 2 and z=0.008 ± 0.004, which 
corresponds to v ou t<0.007c. The fit improvement is Ax 2 =77.8 for three model parameters. This provides the best-fit, 
with F-test probability P F >99.99%. A detection of an Fe K absorber was already reported by Reeves et al. (2004). 
Our results reported in Table 3 are consistent overall with theirs. 

NGC 3516 

Observation 0401210401: the velocity width of the absorption line is a v < 2700 km/s. From Fig. 12 in Appendix 
B, there is only one redshift solution. The initial baseline model statistic is x 2= 263.8 for d.o. f. =160. Using an Xstar 
grid with v t =1000 km/s we obtain log£=4.04i°;^ erg s" 1 cm, N H = (5A± 2 1 °) x 10 22 cm" 2 and z=0.003 ± 0.002, which 
corresponds to v ou t=0.006± 0.002c. The fit improvement is Ax 2 =58.4 for three model parameters. This provides the 
best-fit, with F-test probability P F >99.99%. 

Observation 0401210501: the velocity width of absorption lines is a v < 2800 km/s. From Fig. 12 in Appendix B, 
there is only one redshift solution. The initial baseline model statistic is x 2 =240.3 for d.o. f. =159. Using an Xstar 
grid with w t =1000 km/s we obtain log£=4.18±g;$j erg s" 1 cm, N H =(7.0±l%) x 1q22 cm ~ 2 and ^=0.001 ± 0.002, which 
corresponds to v out =0.008± 0.002c. The fit improvement is Ax 2 =45.3 for three model parameters. This provides the 
best-fit, with F-test probability P F >99.99%. 

Observation 0401210601: the velocity width of the absorption lines is a v < 2600 km/s. From Fig. 12 in Appendix B, 
there is only one redshift solution. The initial baseline model statistic is x 2 =225.1 for d.o. f. =157. Using an Xstar grid 
with u t =1000 km/s we obtain log£=4.05+£^ erg s' 1 cm, N H =(6.7± 2 { 3 4 ) x 10 22 cm" 2 and ^=-0.002 ± 0.002, which 
corresponds to w ou t=0.011 ± 0.002c. The fit improvement is Ax 2 =54.8 for three model parameters. This provides the 
best-fit, with F-test probability P F >99.99%. 

Observation 0401211001: the velocity width of the absorption lines is a v < 2800 km/s. From Fig. 12 in Appendix B, 
there is only one redshift solution. The initial baseline model statistic is x 2 =204.7 for d.o. f =157. Using an Xstar grid 
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with u t =1000 km/s we obtain log£=4.18±^o erg s" 1 cm, N H =(7.ltl\i) x 1022 cnl ~ 2 and ^=-0.004 ± 0.003, which 
corresponds to v O ut=0.013 ± 0.003c. The fit improvement is Ax 2 =25.1 for three model parameters. This provides the 
best-fit, with F-test probability Pp=99.98%. The detection of an Fe K absorber was already reported by Turner et 
al. (2008). Our results reported in Table 3 are consistent overall with theirs. 

Mrk 509 

Observation 0130720101: the velocity width of the absorption line is a v < 3800 km/s. From Fig. 12 in Appendix 
B, there is only one redshift solution. The initial baseline model statistic is x 2 =151.0 for <i.o./.=153. Using an Xstar 
grid with v t =1000 km/s we obtain (1) log£=5.20i£g£ erg s" 1 cm, N H >7 x 10 22 cm" 2 and z=-0.132 ± 0.003, which 
corresponds to v out =0.173± 0.003c. The fit improvement is Ax 2 =H-8 for three model parameters. Using an Xstar 
grid with v t =3000 km/s we obtain (2) log^5.55l?;°? erg s" 1 cm, N H >6 x 10 22 cm" 2 and z=-0.132 ± 0.004, which 
corresponds to f ou t=0.173 ± 0.004c. The fit improvement is Ax 2 =10.6 for three model parameters. The fits (1) and 
(2) are equivalent at 90%. The best-fit is provided by (1), with F-test probability P F =99.4%. We detected a UFO. 

Observation 0306090201: the velocity width of the absorption line is a v < 3200 km/s. From Fig. 12 in Appendix 
B, there is only one redshift solution. The initial baseline model statistic is x 2 =189.7 for c£.o./.=161. Using an Xstar 
grid with u t =1000 km/s we obtain (1) log£=5.73t?;e9 erg s" 1 cm, N H >6 x 10 22 cm" 2 and z=-0.101 ± 0.002, which 
corresponds to u ou4 =0.139 ± 0.002c. The fit improvement is Ax 2 =12.7 for three model parameters. Using an Xstar 
grid with w 4 =3000 km/s we obtain (2) log£=5.73±?;^ erg s" 1 cm, N H >10 x 10 22 cm" 2 and z=-0.100 ± 0.004, which 
corresponds to v OU i=0.138± 0.004c. The fit improvement is Ax 2 =10.1 for three model parameters. The best-fit is 
provided by (1), with F-test probability P F =99%. The fits (1) and (2) are equivalent at 90%. We detected a UFO. 

Observation 0306090401: the velocity width of the absorption line is a v < 5000 km/s. From Fig. 12 in Appendix B, 
there are two possible redshift solutions, corresponding to an identification of the line with Fe XXV Hea or Fe XXVI 
Lya. The initial baseline model statistic is x 2 =215.5 for d.o./.=162. Using an Xstar grid with v t =1000 km/s. For 
Fe XXV we obtain (1) log£=3.47j:o;i5 erg s" 1 cm, N H =(1.0±%i) x 1q22 cm ~ 2 and z=-0.185±0.003, which corresponds 
to v ottt =0.234±0.003c. The fit improvement is Ax 2 =12.5 for three model parameters. Instead, for Fe XXVI we obtain 
(2) log£= 5.13±g;f§ erg s" 1 cm, N H >7 x 10 22 cm" 2 and z=-0.152± 0.003, which corresponds to v out =0. 196 ± 0.003c. 
The fit improvement is A% 2 =14.5 for three model parameters. Using an Xstar grid with f t =3000 km/s. For Fe XXV 
we obtain (3) log£=2.97lg; 2 g erg s" 1 cm, AT H =(0.8±g;f ) x 10 22 cm" 2 and z=-0.188 ± 0.004, which corresponds to 
v out — 0.237 ± 0.004c. The fit improvement is A% 2 =12.2 for three model parameters. Instead, for Fe XXVI we obtain 
(4) log£= 5.54±g;^| erg s" 1 cm, N H >4 x 10 22 cm" 2 and z=-0.153± 0.004, which corresponds to v out =0. 197 ± 0.004c. 
The fit improvement is A% 2 =13.1 for three model parameters. Using an Xstar grid with v t =5000 km/s. For Fe XXV 
we obtain (5) log£=3A8±%^ erg s" 1 cm, N H =(1.0±%%) x 10 22 cm" 2 and z=-0.185 ± 0.007, which corresponds to 
v out — 0.234 ± 0.007c. The fit improvement is Ax 2 =10.1 for three model parameters. Instead, for Fe XXVI we obtain 
(6) log^=5.28^;^ erg s" 1 cm, N H >4 x 10 22 cm" 2 and z=-0. 154 ± 0.005, which corresponds to v out =0. 198 ± 0.005c. 
The fit improvement is Ax 2 =ll-5 for three model parameters. The best-fit is provided by (1), with F-test probability 
Pp=99%. The fits (1), (2), (3) and (4) are equivalent at 90%. The detection of UFOs in these observations was already 
reported by Cappi et al. (2009). Our results reported in Table 3 are consistent overall with theirs. 

Ark 120 

Observation 0147190101: the velocity width of the absorption lines is a v < 3300 km/s. From Fig. 12 in Appendix B, 
there are two possible redshift solutions, corresponding to an identification of the line with Fe XXV Hea or Fe XXVI 
Lyct. The initial baseline model statistic is x 2 =161.8 for rf.o./.=156. Using an Xstar grid with i> t =1000 km/s. For 
Fe XXV we obtain (1) log£=3.44±8;ff erg s" 1 cm, N H =(1.1±%%) x 10 22 cm" 2 and z=-0.247±0.003, which corresponds 
to u otI t=0.306± 0.003c. The fit improvement is Ax 2 =12.5 for three model parameters. Instead, for Fe XXVI we obtain 
(2) log£=5.27±£|? erg s" 1 cm, N H >3 x 10 22 cm" 2 and z=-0. 216 ± 0.003, which corresponds to v out =0.269 ± 0.003c. 
The fit improvement is Ax 2 =11.8 for three model parameters. Using an Xstar grid with f t =3000 km/s. For Fe XXV 
we obtain (3) log£=4.05±?;^ erg s" 1 cm, N H =(1.6±H) x 10 22 cm" 2 and z=-0.247 ± 0.003, which corresponds to 
Wowt=0.306 ± 0.003c. The fit improvement is A% 2 =9.7 for three model parameters. Instead, for Fe XXVI we obtain 
(4) log^=5.68i°;2i er S s_1 cm > N H >4 x 10 22 cm" 2 and z=-0. 216 ± 0.004, which corresponds to v out =0. 269 ± 0.004c. 
The fit improvement is Ax 2 =10.1 for three model parameters. The best-fit is provided by (1), with F-test probability 
Pf=99A%. The fits (1), (2) and (4) are equivalent at 90%. We detected for the first time a UFO in this observation. 

Mrk 279 

Observation 0302480501: the velocity width of the absorption line is a v < 4500 km/s. From Fig. 12 in Appendix 
B, there is only one redshift solution. The initial baseline model statistic is x 2 =146.6 for do./.=155. Using an Xstar 
grid with u t =1000 km/s we obtain (1) log£=4.04±^ erg s" 1 cm, N H =(5.8t.l$) x 1q22 cm ~ 2 and ^=0.031 ± 0.004, 
which corresponds to v out <0.007c. The fit improvement is Ax 2 =14.0 for three model parameters. Using an Xstar grid 
with w t =3000 km/s we obtain (2) log£=4.08t£|j erg s" 1 cm, N H =(3.7±2 5 7 ) x 10 22 cm" 2 and z=0.031 ± 0.005, which 
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corresponds to v otl *<0.009c. The fit improvement is Ax 2 =11.6 for three model parameters. The best-fit is provided 
by (1), with F-test probability P F =99.85%. The fits (1) and (2) are equivalent at 90%. 

Mrk 79 

Observation 0400070201: the velocity width of the absorption line is a v < 6400 km/s. From Fig. 12 in Appendix B, 
there is only one redshift solution. The initial baseline model statistic is x 2 =151.0 for d.o./.=141. Using an Xstar grid 
with v t =1000 km/s we obtain (1) log£=4.17±g;^ erg s^ 1 cm, iV H =(15.2±}. 6 9 2 ) x 10 22 cm" 2 and z=-0.067±0.003, which 
corresponds to v ou t=0.091±0.003c. The fit improvement is Ax 2 =19.1 for three model parameters. Using an Xstar grid 
with v t =3000 km/s we obtain (2) log£=4.34+°;^ erg s" 1 cm, N H =(l%.lt%\) x 10 22 cm" 2 and z=-0.068±0.003, which 
corresponds to u OU i=0.092 ± 0.003c. Ax 2 =16.5 for three model parameters. Using an Xstar grid with y t =5000 km/s 
we obtain (3) log£=4.29t£$ erg s" 1 cm, N H ={ll.btH) x 10 22 cm" 2 and z=-0.068 ± 0.005, which corresponds to 
tw=0.091 ± 0.005c. The fit improvement is Ax 2 =13.5 for three model parameters. The best-fit is provided by (1), 
with F-test probability P F =99.97%. The fits (1) and (2) are equivalent at 90%. We detected for the first time a UFO 
in this observation. 

NGC 4051 

Observation 0109141401: the velocity width of the absorption line is a v < 4200 km/s. From Fig. 12 in Appendix B, 
there are two possible redshift solutions, corresponding to an identification of the line with Fe XXV Hea or Fe XXVI 
Lya. The initial baseline model statistic is x 2 =174.3 for d.o. /.=158. Using an Xstar grid with t> t =1000 km/s. For 
Fe XXV we obtain (1) log£=3.64±g;^ erg s" 1 cm, N H =(1.1±J |) x 10 22 cm" 2 and z=-0.051±0.003, which corresponds 
to v out =0.054 ± 0.003c. The fit improvement is Ax 2 =13.8 for three model parameters. Instead, for Fe XXVI we obtain 
(2) log£=5.73±?;£| erg s^ 1 cm, N H >6 x 10 22 cm" 2 and z=-0.014± 0.005, which corresponds to u out =0.016 ± 0.005c. 
The fit improvement is Ax 2 =ll-8 for three model parameters. Using an Xstar grid with v t =3000 km/s. For Fe XXV 
we obtain (3) log£=4.18±J;§g erg s" 1 cm, N H =(2A±\^) x 10 22 cur 2 and z=-0.054 ± 0.005, which corresponds to 
v ou t— 0-057 ± 0.005c. The fit improvement is Ax 2 =12.0 for three model parameters. Instead, for Fe XXVI we obtain 
(4) log£= 5.69t?;^ erg s" 1 cm, N H >5 x 10 22 cm" 2 and z=-0.015 ± 0.005, which corresponds to w out =0.017± 0.005c. 
The fit improvement is Ax 2 =ll-5 for three model parameters. The best-fit is provided by (1), with F-test probability 
Pf=99.5%. The fits (1), (2), (3) and (4) are equivalent at 90%. Considering the mean value of the outflow velocity 
between the two line identifications in Table 3, this absorber is consistent with a UFO. 

Observation 0157560101: velocity width absorption line is measured as a v = 4700±ifgg km/s. We used an Xs- 
tar grid with u t =5000 km/s. From Fig. 12 in Appendix B there is only one preferred solution. We obtain 
log^=3.06tg;^ erg s" 1 cm, N H =(5.2±\-l) x 10 22 cur 2 and z=-0.183±0.006, which corresponds to u out =0.202±0.006c. 
The initial baseline model statistic is x 2 =174.6 for <i.o./.=134. The fit improvement is Ax 2 =22.0 for three model pa- 
rameters, which corresponds to an F-test probability of Pp=99.95%. We detected for the first time a UFO in this 
observation. 

Mrk 766 

Observation 0304030301: the velocity width of the absorption line is measured as a v = 4600± 2 ™° km/s. We used 
an Xstar grid with y t =5000 km/s. From Fig. 12 in Appendix B there is only one solution. We obtain log£=3. 461^32 
erg s" 1 cm, N H =(2.lt%i) x 10 22 cur 2 and z=-0.067± 0.006, which corresponds to v out =0.082 ± 0.006c. The initial 
baseline model statistic is x 2 =205.0 for d.o. / '.=157. The fit improvement is Ax 2 =22.4 for three model parameters, 
which corresponds to an F-test probability of Pp=99.95%. We detected a UFO. 

Observation 0304030501: the velocity width of the absorption line is a v < 4400 km/s. From Fig. 12 in Appendix B, 
there is only one redshift solution. The initial baseline model statistic is x 2 =208.3 for d.o. /.=154. Using an Xstar grid 
with w t =1000 km/s we obtain (1) log£=4.28t^ erg s" 1 cm, N H ={10.0tsi) x 1q22 cm ~ 2 and ^=-0.072±0.002, which 
corresponds to v ou t—0. 088 ± 0.002c. The fit improvement is Ax 2 =21.2 for three model parameters. Using an Xstar 
grid with w t =3000 km/s we obtain (2) log£=4.40j;^ erg s" 1 cm, N H ={§.2t A 2 %) x 10 22 cm" 2 and z=-0.074 ± 0.003, 

which corresponds to v ou t—0. 090 ± 0.003c. The fit improvement is Ax 2 =17.3 for three model parameters. The best-fit 
is provided by (1), with F-test probability P F =99.9%. The fit (1) is better than (2) at 90% level. We detected for the 
first time a UFO in this observation. The detection of a possible UFO in the first observation was already suggested 
by Miller et al. (2007) and Turner et al. (2007). Our result reported in Table 3 is consistent overall with theirs. 

Mrk 841 

Observation 0205340401: the velocity width of the absorption line is a v < 4000 km/s. From Fig. 12 in Appendix B, 
there are two possible redshift solutions, corresponding to an identification of the line with Fe XXV Hea or Fe XXVI 
Lya. The initial baseline model statistic is x 2 =147.6 for d.o. /.=130. Using an Xstar grid with w t =1000 km/s. For 
Fe XXV we obtain (1) log£=3.46±gJf erg s" 1 cm, N H =(3.5±\-%) x 10 22 cm" 2 and z=-0.035±0.004, which corresponds 
to w ottt =0.071± 0.004c. The fit improvement is Ax 2 =12.2 for three model parameters. Instead, for Fe XXVI we obtain 
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(2) log£=4.50l£^ erg s^ 1 cm, N H >U x 10 22 cm- 2 and z=0.003 ± 0.004, which corresponds to u out =0.034 ± 0.004c. 
The fit improvement is Ax 2 =12.7 for three model parameters. Using an Xstar grid with ^=3000 km/s. For Fe XXV 
we obtain (3) log£=2.93±o;°7 erg s" 1 cm, N H =(2.3±f%) x 10 22 cm" 2 and z=-0.038 ± 0.005, which corresponds to 
v out ~0.075 ± 0.005c. The fit improvement is Ax 2 =11.9 for three model parameters. Instead, for Fe XXVI we obtain 
(4) log£=4.54±°;°^ erg s" 1 cm, N H >10 x 10 22 cm" 2 and z=0.003 ± 0.005, which corresponds to *w=0.034 ± 0.005c. 
The fit improvement is A% 2 =12.4 for three model parameters. The best-fit is provided by (2), with F-test probability 
Pf=99.07%. The fits (1), (2), (3) and (4) are equivalent at 90%. We detected for the first time a UFO in this 
observation. 

ESO 323- G77 

Observation 0300240501: the velocity width of the absorption lines is a v < 3700 km/s. From Fig. 12 in Appendix B, 
there is only one redshift solution. The initial baseline model statistic is x 2 =223.8 for d.o. /.=135. Using an Xstar grid 
with w t =1000 km/s we obtain (1) log£=3.45t£ra erg s" 1 cm, N H =(13.3±H) x 10 22 cm- 2 ' and z=0.004± 0.003, which 
corresponds to u o «t=0.011 ± 0.003c. The fit improvement is Ax 2 =74.3 for three model parameters. Using an Xstar 
grid with v t =3000 km/s we obtain (2) log£=4.03i£^ erg s" 1 cm, N H =(29.1±f 9 2 ) x 10 22 cm" 2 and z=0.003 ± 0.003, 
which corresponds to v out ~0. 012 ± 0.003c. The fit improvement is Ax 2 =77.0 for three model parameters. The best-fit 
is provided by (2), with F-test probability Pf>99.99%. The fits (1) and (2) are equivalent at 90%. The detection of 
an Fe K absorber was already reported by Jimenez-Bailon et al. (2008). Our results reported in Table 3 are consistent 
overall with theirs. 

1H0419-577 

Observation 0148000201: the velocity width of the absorption line is a v < 7000 km/s. From Fig. 12 in Appendix B, 
there is only one redshift solution. The initial baseline model statistic is x 2 =128.7 for do./.=105. Using an Xstar grid 
with w t =1000 km/s we obtain (1) log£=3.85±o ^ erg s" 1 cm, N H =(17.2tl% 3 ) x 10 22 cm" 2 and z=0.023± 0.004, which 
corresponds to i> out =0.076 ± 0.004c. The fit improvement is Ax 2 =14.8 for three model parameters. Using an Xstar grid 
with w t =3000 km/s we obtain (2) log£=3.89±^ erg s" 1 cm, A^=(11.3tg° 3 7 ) x 10 22 cm~ 2 and z=0.019± 0.005, which 
corresponds to v ou t=0.080 ± 0.005c. The fit improvement is Ax 2 =12.2 for three model parameters. Using an Xstar 
grid with u t =5000 km/s we obtain (3) log£=3.70tg;^ erg s" 1 cm, N H =(4.2tH) x 10 22 cur 2 and z=0.021 ± 0.006, 
which corresponds to v out ~0- 078 ± 0.006c. The fit improvement is Ax 2 =ll.l for three model parameters. The best-fit 
is provided by (1), with F-test probability Pp>99.43%. The fits (1) and (2) are equivalent at 90%. We detected for 
the first time a UFO in this observation. 

Mrk 290 

Observation 0400360601: the velocity width of the absorption line is a v < 4000 km/s. From Fig. 12 in Appendix B, 
there are two possible redshift solutions, corresponding to an identification of the line with Fe XXV Hea or Fe XXVI 
Lya. The initial baseline model statistic is x 2 =103.7 for d.o. /.=104. Using an Xstar grid with v t =1000 km/s we obtain 

(1) for Fe XXV log£=3. 54^^25 erg s" 1 cm, N H ={^-\i) x 1q22 cm ~ 2 and z=-0. 142 ± 0.004, which corresponds to 
v o «i=0.181 ± 0.004c. The fit improvement is Ax 2 =11.6 for three model parameters. Instead, for Fe XXVI we obtain 

(2) log^=4.37tg;g erg s" 1 cm, N H >ll.b x 10 22 cm" 2 and z=-0.107±0.004, which corresponds to w out =0.142±0.004c. 
The fit improvement is Ax 2 =13.3 for three model parameters. Using an Xstar grid with y t =3000 km/s we obtain 

(3) for Fe XXV logC=3.07ig;^ erg s" 1 cm, V H =(3.lt}; 2 ) x 10 22 cur 2 and z=-0. 143 ± 0.005, which corresponds 
to w otl t=0.182 ± 0.005c. The fit improvement is Ax 2 =12.0 for three model parameters. Instead for Fe XXVI we 
obtain (4) log£=4.35t^ erg s" 1 cm, AT H =(33.3in^) x 10 22 cm" 2 and ^=-0.107 ± 0.004, which corresponds to 
^W=0.142 ± 0.004c. The fit improvement is Ax 2 =13.5 for three model parameters. The best-fit is provided by (4), 
with F-test probability P F >99.74%. The fits (1), (2), (3) and (4) are equivalent at 90%. We detected for the first time 
a UFO in this observation. 

Mrk 205 

Observation 0124110101: the velocity width of the absorption line is a v < 3100 km/s. From Fig. 12 in Appendix 
B, there is only one redshift solution. The initial baseline model statistic is x 2 =135.7 for d.o./.=121. Using an Xstar 
grid with u t =1000 km/s we obtain (1) log£=4.86to;go er S s_1 cm : N H >U.5 x 10 22 cur 2 and z=-0.031 ± 0.004, which 
corresponds to iw=0.100 ± 0.004c. The fit improvement is Ax 2 =H-9 for three model parameters. Using an Xstar 
grid with u 4 =3000 km/s we obtain (2) log£=5.3lt^ erg s" 1 cm, N H >22.7 x 10 22 cur 2 and z=-0.031 ± 0.004, which 
corresponds to v ou t =0. 100 ± 0.004. The fit improvement is Ax 2 =10.1 for three model parameters. The fits (1) and (2) 
are equivalent at 90%. The best-fit is provided by (1), with F-test probability Pp>98.8%. We detected for the first 
time a UFO in this observation. 
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PG 1211+143 

Observation 0112610101: the velocity width of the absorption line is measured as a v = 4500+*2oo km/s. We used an 
Xstar grid with v t —5000 km/s. From Fig. 12 in Appendix B, there is only one solution. We obtained log£=2.87tQ':[ 2 1 
erg s" 1 cm, N H =(8.0±f%) x 10 22 cm" 2 and z=-0.071 ± 0.003, which corresponds to v out =0. 151 ± 0.003c. The initial 
baseline model statistic is x 2 =196.0 for d.o.f. =117. The fit improvement is Ax 2 =70.5 for three model parameters, 
which corresponds to an F-test probability of Pp>99.99%. We detected a UFO. The detection of a UFO in this 
observation was already reported by Pounds et al. (2003) and Pounds & Page (2006). Our results reported in Table 3 
are consistent overall with theirs. 

MGG-5-23-16 

Observation 0302850201: the velocity width of the absorption line is a v < 13000 km/s. From Fig. 12 in Appendix B, 
there is only one redshift solution. The initial baseline model statistic is x 2 =298.4 for d.o.f .=161. Using an Xstar grid 
with v t =1000 km/s we obtain (1) log£=4.30±g;?g erg s" 1 cm, N H =(3.2±%%) x 10 22 cm" 2 and z=-0. 101 ± 0.003, which 
corresponds to iw =0.1 15 ± 0.003c. The fit improvement is Ax 2 =34.0 for three model parameters. Using an Xstar grid 
with v t =3000 km/s we obtain (2) log£=4.38±g;$j erg s" 1 cm, A#=(4.3tj]-j) x 10 22 cm" 2 and z=-0.103 ± 0.003 which 
corresponds to v o „i=0.117± 0.003c. The fit improvement is A% 2 =41.3 for three model parameters. Using an Xstar 
grid with v t =5000 km/s we obtain (3) log£=4.31±[};$j erg s" 1 cm, A^=(3.8t}^) x 10 22 cm" 2 and z=-0.102 ± 0.004, 
which corresponds to v out —0. 116 ± 0.004c. The fit improvement is Ax 2 =41.0 for three model parameters. The best-fit 
is provided by (2), with F-test probability Pp>99.99%. The fits (2) and (3) are equivalent at 90%. The detection of a 
UFO in this observation was already suggested by Braito et al. (2007). Our results reported in Table 3 are consistent 
overall with theirs. 

NGC 4507 

Observation 0006220201: the velocity width of the absorption line is a v < 4200 km/s. From Fig. 12 in Appendix 
B, there are two possible redshift solutions, corresponding to an identification of the line with Fe XXV Hea or 
Fe XXVI Lya. Using an Xstar grid with v t =1000 km/s we obtain (1) for Fc XXV log£=3.81±°;|3 erg s" 1 cm, 
N H = (2.3t 2 1 8 3 ) x 10 22 cm" 2 and z=-0.188 ± 0.003, which corresponds to v out =0. 216 ± 0.003c. The fit improvement 
is Ax 2 =10.1 for three model parameters. Instead, for Fe XXVI we obtain (2) log^=5.36^Q;g 2 erg s" 1 cm, A#>10 x 
10 22 cm and z= — 0.155 i 0.003, which corresponds to v ou i — 0.178 ± 0.003c. The fit improvement is Ax 2 =8.7 for 
three model parameters. Using an Xstar grid with v t =3000 km/s we obtain (3) for Fc XXV log£=4.16+o;7g erg s" 1 cm, 
N H =(4.1±f%) x 10 22 cm" 2 and z=-0.190 ± 0.004, which corresponds to w otlt =0.218 ± 0.004c. The fit improvement is 
Ax 2 =9.7 for three model parameters. Instead, for Fe XXVI we obtain (4) log£=4.78to;3y erg s" 1 cm, Nh >5x 10 22 cm" 2 
and z=— 0.156 ± 0.004, which corresponds to w ott t=0.179 ± 0.004c. The fit improvement is Ax 2 =8.3 for three model 
parameters. The best-fit is provided by (1), with F-test probability P F =97.7%. The fits (1), (2), (3) and (4) are 
equivalent at 90%. We detected for the first time a UFO in this observation. 

NGC 7582 

Observation 0112310201: the velocity width of the absorption line is a v < 11000 km/s. From Fig. 12 in Appendix B, 
there is only one redshift solution. The initial baseline model statistic is x 2 =131.9 for d.o./.=103. Using an Xstar grid 
with v t =1000 km/s we obtain (1) log£=3.39±£?§ erg s" 1 cm, N H =(2SAt 2 9 1 s 1 ) x 10 22 cm" 2 and z=-0.250±0.003, which 
corresponds to i> O ut=0.285±0.003c. The fit improvement is Ax 2 =20.9 for three model parameters. Using an Xstar grid 
with v t =3000 km/s we obtain (2) log£=3.74tg;jg erg s" 1 cm, N H =(U.0+l% 2 ) x 10 22 cm" 2 and z=-0.250±0.003, which 
corresponds to w OU i=0.285 ± 0.003c. The fit improvement is Ax 2 =18.1 for three model parameters. Using an Xstar 
grid with v t =5000 km/s we obtain (3) log£=3.31±o j§ erg s" 1 cm, N H =(6.&tH) x 1q22 cm ~ 2 and ^=-0.248 ± 0.005, 
which corresponds to v out =0. 283 ± 0.005c. The fit improvement is Ax 2 =15.6 for three model parameters. The best-fit 
is provided by (1), with F-test probability P F =99.94%. The fit (1) is better than (2) and (3) at the 90% level. We 
detected for the first time a UFO in this observation. 

We note that TA reported the simultaneous evidence for an absorption line at E~9 keV, possibly due to Fe XXV, and 
three additional weak lines in the energy range E~4-5.3 keV, possibly due to Hc/H-likc Ar and Ca, in this observation. 
The significance of the lower energy lines is only marginal, being in the range ~95-99%, and their modeling can be 
hampered by the intrinsic complexity of the spectrum of this source at energies lower than 6 keV, showing the 
superimposition of different emission/absorption components (Piconcelli et al. 2007). However, as already discussed in 
TA, if the centroids of the lower energy lines are fixed to the expected energies for Ar XVII Hea (E~3.1 keV), Ar XVII 
Lya (E~3.3 keV) and Ca XIX Hea (E~3.9 keV) (Verner et al. 1996) and they are fitted together with the Fe XXV 
Hea line, it turns out that they are consistent with sharing the same common blue-shift of z =—0.255 ± 0.002. This 
modeling provides a good representation of the data, with a fit improvement of Ax 2 =35 for five parameters, which 
corresponds to a very low random probability of less than 10" 5 . This was already discussed in §4.5 and Appendix B of 
TA. However, from the ratios of the EWs of the lower energy lines with respect to the Fe XXV Hea we can also roughly 
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estimate that there is the need for a ^5 overabundance of lighter elements, such as Ar and Ca, with respect to Fe. In 
particular, even though in this case the statistics is presumably mostly driven by the intense Fe XXV/XXVI absorption 
lines, a single Xstar grid seems capable to simultaneously parametrize all the lines, but those at lower energy are more 
intense than what expected assuming standard Solar abundances. An Xstar grid with ~5 overabundance of Ar and 
Ca actually provides a better modeling of these lines. However, the detailed analysis of possible low energy lines is 
beyond the scope of the present paper, which is focused on the study of the UFOs essentially through the more intense 
Fe K-shell lines. We postpone their treatment to forthcoming works. 
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X 2 STATISTIC AND OBSERVED ABSORBER REDSHIFT CONTOUR PLOTS 
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Figure 12. \ statistic plotted against the Xstar absorber redshift. The black solid, dashed red and dotted blue contours refer 
to Xstar grids with turbulent velocities of 1000 km/s, 3000 km/s and 5000 km/s, respectively. Not all observations required 
grids with different turbulent velocities, see text for more details. The horizontal lines indicate the 90% confidence levels for one 
interesting parameter. The vertical dash-dotted line indicates the cosmological redshift of the source. The remaining vertical 
lines indicate the absorber redshift from the Xstar fits. In this case a negative redshift indicates a blueshift. 
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Figure 14. Continued. 



